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ABSTRACT 


This manual describes how to use the Emulation/Simulation Computer Model, 
ESCM Based on G189A, .ESCM computes the transient performance of a Space 
Station atSpUe revitalization subsystem (ARS) with C0 2 removal provided 
by a solirS^water desorbed subsystem called SAWD Many performance 
parameters are computed some of which are cabin C0 2 partial pressure, relative 
humicM ty^temperatu re , 0 2 partial pressure, and dew point. The program allows 
the user to simulate various possible combinations of man loading, meta o 
profits, cabin volumes and certain hypothesized failures that could occur. 
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1.0 

INTRODUCTION 

l.i 

Background 


1.2 


„ whirh was developed to demonstrate the 

C Xr po^on S’fhe^uuop/Ulat^on Si ting, and 
uwjui ^ /ccctsd^ See Reference 1. ESbiAH is 


ESCM is a 
utility of a iiwjut 

Technology Assessment Program (ESS TAP ^“Xite'm.rengineerlng 

a concept for software tools that will support tner a 

process" of the Space Station begihnin; I with T «"«Pj e ^ tion of the 

continuing on through mission op nnora ti 0 nal chases of a flight 
software tools with the design and < .per *^1 Ph*se hl ° 5 sof tware 

system is shown in Figure 1. The P h ^“ s W 0 1n the design 

concept is to have the ord er to provide the 


process precede development i 
greatest design and cost benefits. 

cnaro station the Environmental Control 
Of the many systems in theS P a selected to demonstrate the 

skis, •«. : *?r ‘.as 

program was targeted to eva1 ^ e a " d ^posea" evaluation of the 
phases 2 through 4 software. For this purpose, an ev program 

T "r~h C cT* f T'°Z “a a re?uuTthis thj 

Su-bs,stem 

desi gn ^tool "ut^fuy^and the independent nature of the ARS compared 
to other subsystems. 

The resulting ESCM model has been ^ ^ 

computer program in the design, ev P .’ the use of t he program, 

development hardware. Specifically, ified Failure Mode and 

design hardware and software planning can he 

Effect Analysis (FMEA) can De asbi^eu, 

technically improved. 

Computer Program Overview 

The ESCM Computer Pr ?>S r ** ^ s J*** s ?e r f 'progra^ EOT« thellSSA 
3. Like a generalized heat transfer p ^ G189A subroutine 

library"" for "model's .1, o/the components and 1 lastly, the 

G189A executive routines to take the input u the 

SST?r™ ■* "the S GP0LY1 and” GPOLyV Subroutines which 
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Figure 1 

Application Of Software Tools To The Design 
And Operational Phases Of A Flight System 
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Figure 2 

ARS - SAWD II Schematic For Demonstration Model 
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Figure 3 

Overview Of ESCM Computer Program 
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2.0 


Computer Program Overview (Continued) 


define the control logic for. the 

developed component ^ ubr0 ^ tin ® . hich nrovides computer generated 
'r^t^^o^lhf^Lh^^p^^th^^orJanrSrrfo^nce palters 
Ind ubular printoff to a terminal “ r “" f ® r 'subgroup for 

G 189 A^ 4 and^^hich* 1 * contains 1 ^many^ 1 performance constants and inpot 
parameters . 


ESCM starts with input data which J®f’ ne ? t the w J5J 
ents and plumbing and combines this data 
options. Typical input parameters are: 


(1) Number of men 

(2) Cabin volume 

(3) Cabin temperature 


(4) Equipment heat 

(5) Cabin leakage 

(6) Metabolic profiles 


subgroup compon- 
various program 


load 


ESCH then solves each of the & VS^E". 

Hr SMlS* X? typical perfo^ nee parsers 
to the screen and a schematic £ ^/^/“^/the frequency of 
printoff . P 'lypCca 1* computed performance parameters are: 

(1) Cabin CO 2 mmHg 

(2) Cabin O 2 Partial pressure 

(3) Cabin relative humidity 

(4) Cabin temperature 

(5) Cabin total pressure 

After the display to the screen, the ^ rime step 

GP0LY1 is solved, all th % c Q % P 2 ° ne s solved and lastly ARS subroutine 
the specified sequence GP0LY2 1 d ormance parameters. 

The process 'is ^ repeated 3 "unt^V^ he ^time specified by the user is 

reached . 

description OF HARDWARE AND MODEL 

ESCM is a computer p rogram , whi ch - J™ 1 ^ECLSS f.^To demon- 

Envi ronmental Control and Li e PP FCLSS had to be modeled. 

sm, s Ts iy a 

LeliranT t 9 i°ve W a''de«°Hp d tio a n “"the ESCM computer program. 
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2.1 Description of System to be Modeled 

As stated above, the system modeled is a subgroup of the 
realization subsystem. The major elements are the cabin, the 
crew V sensible hlat exchanger, a humidity control "eat exchanger, 
fans, a SAWD carbon dioxide removal system, and a two g 

trol ler. 

The cabin is a volume of user specified size which houses the crew, 
two gas controller, and certain user specified vehicle loads 
representative of test equipment. The user may also specify air 
leakage out of the cabin. 

tl 0 -raw nf user specified size and metabolic loading breathes in 
oxygen 6 f ran he ca'bln’ air and releases carbon dioxide an wa er 
vapor. The crew also gives off sensible heat to the cabin air. 


heat 

heat 


Cabin air is drawn off by one fan to pass th e aj; through^ 

exchanger “Ehfa™ heat* to™, ant wat^r which enters 

Sp^^1^:,Ur^in r coKt^I 1 ?lie t ^n¥^te ^‘maintain 
the desired cabin air temperature. The heat exchanger p 

fin type similar to that on the Shuttle. 

™ ss.“nrr^sJ5? sr. S-rrr; 

^^this heat exchanger. A portion of the mam air stream with 
condensed vapor is sent to a water separator where the c ° nde ^® 
“p^r is !^- d i 

system^ Ito'lE SfSV.“ 2? store from the SAWO system on 
cabi n reldtive humidity is lessened* 

TL. last f an draws air from the cabin for use by the SAWD carbon 
Hinvi do removal system. This carbon dioxide removal system contains 
?wo Imlne K Solenoid valves, an accumulator, a water pump steam 
generator, a controller, and necessary plumbing and 1 ^^"co^tin- 

LrS ’ alternately 11 a& *UT^" ‘ "™e 

zrkkrsir « tzx s , «M. rt Sd e s 1 n » 

s?andby Initial flow from the desorbing bed 12 flows to the cabin 
until Y the flow reaches 0.01 cfm as detected by the Tlow sensor. 
Then the flow is directed to an accumulator. Pressure builds in 
the accumulator until 24 psia is reached. At this pressure flow 
nprmitted to leave at a controlled rate to a carbon dioxme 
reduction process. If pressure builds beyond 30 psia, carbon 
dioxide is dumped overboard to lower the f res ^ 

During desorb, the pressure in the accumulator back pre: s 
desorbing bed. A check valve prevents any back flow from t 
accumulator to the beds. 
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2.1 


Description of System to be Modeled (Continued) 

Desorbinq in bed #2 continues until its exit temperature reaches 

temperature reaches 130°F, whichever occurs first. 

After energy transfer is completed, steam begins to flow to bed #1 

Llrb Tt while valve positions are changed to direct air flow 
to desorb it wniie vaive p Desorbing continues in bed #1 

Sms!?’ ^ 

*JS» s?. 

sequence as directed by the SAWD controller. 

Air from the SAWD system enters the ^mldUy control heat exchanger 
where any moisture is removed. The air returning 
the SAWD system now has less carbon dioxide. 

Cabin air total pressure and oxygen partial pressure are controlled 
Shuttle” 0 

to£ pressure H i a and the oxygen partial pressure at 3.09 
)n 3 23 osia Details of the control logic are shown in Table 2. 
The Opening and closing curves are input by the user as part of the 
input table data. See Section 4.0. 

2.2 Description of Model 

The previously described system is modeled for - use in the ESCM 

ra :• ^ *3 

paths °aTe °gfven a $ letter" P* o^fto’de^gnate tSe^ath into or out of 
the component as primary or secondary. 

Tahlp 3 identifies the component subroutine used for each component 

and also presents the G189A subroutine number All components are 

labelled in Figure 4 except t «. enuring ani one 

r. ngTthe com, onent is a mixer. It simply mixes the two entering 
flows. 9 If the small circle has one entering flow and wo eaving, 

the component is a splitter. It splits tne ,, the 

two paths per a user provided split ratio. Of course, all 
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Table 2 


Controller Logic For H.5 psi 


Shuttle Two Gas Controller 


Initial 

Oxygen 

Valve 

Status 

N/A 

N/A 

Closed 

Open 

Closed 

Open 

N/A 

N/A 

N/A 

N/A 


-Condition 


■Action 


Initial 

Nitrogen 

Valve 

Status 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

Closed 

Open 

Closed 

Open 


Last 

Valve 

Opened 

PT 

( psi a ) 

PO 2 

(psia) 



N/A 

>14.819 

— 

N/A 

>14.813 

>3.23 

N/A 

<14.819 

<3.09 

N/A 

<14.819 

<3.09 

Oxygen 

<14.819 

3.09<P02<3.23 

Oxygen 

<14.819 

3.09<P02<3.23 

Nitrogen 

<14.819 

3.09<P02<3.23 

Nitrogen 

<14.819 

3.09<P02<3.23 

N/A 

<14.813 

>3.23 

N/A 

<14.813 

>3.23 


Oxygen 

Valve 

Status 

Nitrogen 

Valve 

Status 

0 2 Flow Per 
Opening Or 
Closing Curve 

Closed 

Closed 

— 

Closed 

Closed 

— 

Open 

Closed 

Opening 

Open 1 

Closed 

Closing 1 

Open 

Closed 

Openi ng 

Open 1 

Closed 

Closing 1 

Closed 

Open 

— 

Closed 

Open 1 

— 

Closed 

Open 


Closed 

Open 1 

— 


1 When 


flow calculated by closing curve 


9503 


N 2 Flow Per 
Opening Or 
Closing Curve 


Opening 

Closing 1 

Opening 

Closing 1 


= 0.0, valve closes. 
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Figure 4 

ARS - SAWD II Schematic As Modeled In The ESCM 
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ESCM Component Number 


Table 3 

Subroutine Cross Reference Table 


Component 

No, 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


Component Desc ription 


;rew in cabin 
:abin 

Split to humidity control HX or 5AWU 
Cabin Sensible HX fan 
Cabin sensible HX 

Cabin sensible HX cooling fluid boundary 
Cabin condensing HX 
Condensing HX cooling fluid boundary 
Split to water separator or bypass 
Water separator 

Mix water separator and bypass 

Condensing HX fan 

Mix SAWO and condensing HX flows 

Split to SAWO bed #1 or #2 

Mix steam or cabin air to bed #1 

SAWO bed #1 

Split to cond. HX or C02 recovery 
Mix bed #1 and bed #2 flows 
SAWO fan 

Mix steam or cabin air to bed #2 
SAWO bed #2 

Split to cond. HX or C02 recovery 
Mix bed #1 and bed #2 exit gases 
Split to preheat or C02 accumulator 
Mix steam generator and preheat streams 
Split to SAWO bed #1 or #2 
Water supply tank 
Water supply pump 
Steam generator 

Split desorbed gas to cabin or C02 accumulator 
Mix cabin air and SAWO air 
Carbon dioxide accumulator tank 


Subroutine 
Name 


SUITS 

CABIN 

SPLIT 

FAN 

ANYHX 

ALTCOM 

ANYHX 

ALTCOM 

SPLIT 

ALTCOM 

GASMIX 

FAN 

GASMIX 

SPLIT 

GASMIX 

IR45 

SPLIT 

GASMIX 

FAN 

GASMIX 

IR45 

SPLIT 

GASMIX 

SPLIT 

GASMIX 

SPLIT 

TANKG 

PUMP 

SMGEN 

SPLIT 

GASMIX 

TANKG 


Subrouti ne 
No. 


2 

1 

10 

23 

4 

49 

4 

49 

10 

49 

6 

23 

6 

10 

6 

73 

10 

6 

23 

6 

73 

10 

6 

10 

6 

10 

30 

22 

27 

10 

6 

30 


n 
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2.2 Description of Model (Continued) 

entering flow could be directed to take one of the exiting flow 
paths. Splitters and mixers are used to represent valve functions. 
The small square components are alternate components and simply 
provide boundary condition water flows and temperatures to the two 
heat exchangers. 

3.0 DESCRIPTION OF ESCM COMPUTER PROGRAM 

The total ESCM program consists of eight entities: 

• G189A executive subroutines 

• G189A utility subroutines 

• G189A special purpose subroutines 

• G189A plot subroutines 

• G189A component subroutines 

• G189A application specific subroutines GP0LY1 and GP0LY2 

• G189A subroutine additions for ESCM: IR45, ARSS 

• G189A input specific for ESCM 

While a detailed discussion of the sequence of operation G189A can 
be found in Reference 3 for the interested user, a brief description 
is presented here and shown in Figure 5. 

(1) At the start of execution, the G189 executive subroutines 
process the input data by allocating the data storage space and 
dynamically loading the integer and floating point data into a 
large single array; the K and V array. 

(2) After the input data is processed, a printoff of the component 
initial values and the solution path is generated. 

(3) From the input solution path, the executive subroutines 
retrieve and unpack from the large single array the first 
component's integer instruction data; and the executive 
subroutines retrieve from the large single array and store into 
small working arrays the component floating point data. 

(4) Thermophysical property data (specific heat, molecular weight, 
thermal conductivity, and viscosity) for the component source 
streams are evaluated using G189A utility subroutines. 

(5) A transfer to the user coded subroutine GP0LY1 is made to allow 
data modifications or logic changes to be incorporated prior to 
component solution. 

(6) The executive routines then transfer computer operation to the 
proper subroutine for the component. 

(7) Following component solution, a transfer to the user coded 
subroutine GP0LY2 is made to allow data modifications or logic 
charges after component solution. 
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Figure 5 

Flow Chart Of ESCM Computer Program Operation 
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(8) The new component status and data are returned to the large 
single array. 

(9) Steps (3) through (8) are repeated for each component in the 
solution path. 

(10) A printoff of the ARS subgroup schematic is made by subroutine 
ARSS for every steady state pass and after every user spec- 
ified number of time steps in transient runs. 

(11) Steps (3) through (10) are repeated until a steady state 
solution is achieved; then, if specified, a transient begins 
and will end at the user specified time. 


The following sections present a description of the seven entities 
previously mentioned. 


3.1 


G189A Executive Subroutines 


The executive subroutines consist of the master control block 
subroutine ECLST and the input editor subroutines: IEDIT, CASDAT, 
BASDAT, MERGEC, KVREAD, TABLRD, SCRUP, ENKODE/DEKODE , IVBLOK, ULOSE , 
UWIN , and DBLOCK. 


The main program is called X189#A and simply transfers operation to 
the master control block subroutine ECLST. 


Complete descriptions of these subroutines can be found in Refer- 
ence 3. 


3.2 G189A Utility Subroutines 

Utility subroutines are those subroutines which are general enough 
in form and application that they could be used by a number of 
subroutines. The utility subroutines used by ESCM are presented in 
Table 4; a complete description of these subroutines can be found in 
Reference 3. 

3.3 G189A Special Purpose Subroutines 

Special purpose subroutines are lower level subroutines which are 
generally called by either the executive subroutines or by a compon- 
ent subroutine. These subroutines were created for one particular 
purpose and are not general enough to be classified as a Utility 
Subroutine. The special purpose routines used by ESCM are presented 
in Table 5; a complete description of these subroutines can be found 
in Reference 3. 

3.4 G189A Component Subroutines 

Each particular type of component in an Environmental Control 
and Life Support System is simulated by a component subroutine. 
The G189A component subroutines used by ESCM are presented in Table 
6, a complete description of these subroutines can be found in 
Reference 3. 
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Table 4 

G189A Utility Subroutines 


Subroutine 

Name 

abrprt 

CYCLE 

ELAPSE 

ESTIM 

FL0ARY 

FL(#SUM 

HBALNC 

HF 

HG 

KK 

LK 

LV 

PR0P 

PRTIME 

PSAT 

QSURR 

RH 

SHELL 



, B, and R Array Printout 
ycle Displacement Time 

; ent ral Processor and Peripheral Processor Elapsed Times 

^estimator for Iterative Calculations 

Selects and Stores Component Flow Data into 19 Member 

Sums Constituent Plows to Obtain Total Flows within a 19 Member 
Flow Array 

* iq Mpmber Flow Array 

Computes Adiabatic Mixture Conditions 
Computes Specific Enthalpy of Liquid Water 
Computes Specific Enthalpy of Water Vapor 

Locates and Returns Value Stored in Specified Component K Array 
Reference Location 

Calculates and Returns K Array Address of Specified Component 
Array Reference Location 

Calculates and Returns V Array Address of Specified Component 
Array Reference Location 

Calculates Thermo-Physical Property Data for a 19 Member Flow 
Array 

Accesses and Prints Computer Clock Time 

Determines Water Vapor Pressure at a Specified Temperature 
Determines Heat Loss to the Surroundings for a Component 
calculates Relative Humidity and Dew Point of a 19 Member Flow 
Array 


npacks Modular K Array Data 
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Table 4 (Continued) 
G189A Utility Subroutines 


Subroutine 

None 


Description 


SK 

ST0PIT 

SV 

TAPE IT 
TSAT 

VALUE 

VARPRT 

VDATPT 

VV 

WASP 


Stores Integer Data into a Specified Component K Array Reference 
Location 

Places End-of-File Mark on Output Tape (Normal Run Termination 
Subroutine) 

Stores Floating Point Data into a Specified Component V Array 
Reference Location 

Used to Dump D/V Array Data onto File 

Determine Saturation Temperature of Water Vapor at a Specified 
Pressure 

Table Interpolation Function 

Labelled Printout of all Components' K and V Array Data 
Component's V Array Data Printout 

Locates and Returns Value Stored in Specified Component V Array 
Reference Location 

Calculates Thermodynamic and Transport Properties of Water and 
Steam 
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Table 5 

G189A Special Purpose Subroutines 

Subroutine 


None 

Description 

ALFRED 

Error Mode Output Format 

C0MS0L 

Component Subroutine Selector 


C02CP 

Computes Effective Specific Heat for Precipitation of Carbon 
Dioxide in a Heat Exchanger 

EFFCP 

Effective Specific Heat for Heat Exchanger 

Cal culations 

HXPER 

Steady State Heat Exchanger Performance 


VLH20 

Calculates Equilibrium Pressure of Water 
Amine Resin 

Vapor on I R -45 Solid 
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Table 6 

G189A Component Subroutines 


Subroutine 

None 

Description 

ALTC0M 

Alternate Component 

ANYHX 

Compact Heat Exchanger 

CABIN 

Cabin or Compartment 

DUCT 

Length of Gas Circuit Ducting 

FAN 

Fan, Blower, or Compressor 

GASMIX 

Gas Mix 

LIQMIX 

Liquid Mix 

PUMP 

Pump 

SMGEN 

Steam Generator 

SPLIT 

Generalized Split 

SUITS 

Crewmen in Suits or Cabin 

TANKG 

Storage Tank 

VACPMP 

Vacuum Pump 
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3.5 


3.6 


G189A P lot Subroutines 

Hat* for use by a separate plotting 
^ t f„ b ™ U V^\ r u\™tines S :?e MAINPL and PLEDIT both of which am 

described in Reference 3. 

GP0LY1 and GP0LY2 for ESCM 

GP0LY1 and SP01Y2 .are coded for ^^‘^^/p^'chanses 1 . 0 GPOLYMs 
executed*)^/ ore °the "^otine while GP0LY2 erecoted 

after the component subroutine. 

Both GPOLY subroutines are setup in the same manner, 
construction is: 


The standard 


TC /M fjp 


n\ r_n T(1 3QQ 


399 CONTINUE 

N is the component number which is about • to ^r just was \" a ^ 

This construction causes logic to be P ti0 „ of this type is 

written ^for^ach component^here desi red by the user. 

3.6.1 GP0LY1 essentially simulates the control 

^ctionVVe"’ror ln the APS subgroup modelled. These control 
functions are: . . . 

(1) ^eH of oZfn partia^pressur^and^cabin”?^^ pressure. 

(2) control of cooling flow to sensible heat exchanger to maintain 
desired cabin temperature. 

,3) control cycling and total operation of SAWD carbon dioxide 
removal system. execut1on of GP0LY1: 

The following are input once at 


DTI ME 
VV( 2,184) 
VV( 16 ,93) 
VV(16 ,94) 
VV( 16 ,79) 
VV( 16 ,84) 
VV( 16,78) 
VV(16,66) 


PHntoff P frequency, s Jj ep # 5 P ib m Print ° ff 

Spec! f *ic^heat°of f foam' in bed « Btu/lbm-F 

: ^c!fn n inro) h ^;l-ih;^« 

; ’r«in in bed #1. Btu/lbm-F 
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3.6.1 


GP0LY1 (Continued) 

The following are input and output from GP0LY1 before solution of 
component number 1 (the crew): 

Input : 

TIME = Time in transient, seconds 

VV(2,104) = Cabin temperature, °F 


Output : 

R(66) = Crew total sensible heat, Btu/hr 

R ( 67) = Crew total latent heat, Btu/hr 

R(82) = Crew total metabolic rate, Btu/hr 

The following are input and output from GP0LY1 before the solution 
of component number 2 (the cabin): 


Input : 
STEADY 
KSYSPAS 


I FREQ 

DTIME 

R(4) 

R(94) 
R( 180) 


= Logical variable denoting if steady state calcula- 
tions are now being performed or not. 

= Integer variable which counts the number of passes 
of the computer program through the solution of all 
the components. 

= The number of computer passes to be performed per 
printoff using ARSS subroutine. 

= Current time step, seconds. 

= Total pressure of flow leaving cabin via the prim- 
ary flow path. 

= Cabin partial pressure of oxygen, psia. 

= Flag denoting whether oxygen or nitrogen was last 
gas admitted through two gas controller. 1 = 
nitrogen, 0 = oxygen. 


Output : 

R ( 165 ) = Flow oxygen into cabin, lb/hr 

R ( 166 ) = Flow of nitrogen into cabin, lb/hr 

The following are input to and output from GP0LY1 before the solu- 
tion of component number 3 (Splitter for flow from cabin): 
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3.6.1 GP0LY1 (Continued) 


Input : 
STEADY 
KSYPAS 

VV( 19 ,1 ) 

V V ( 24 , 1 ) 

VV( 12,1) 

VV( 10,67) 

Output : 
R(65) 


= Logical variable denoting if steady state calcu- 
. number of passes 

of the computer program through the so u to 

- Total °qas n flow ‘through SAWD fan, lb/hr. 

I Total ?low to splitter from SAWD beds, lb/hr. Flow 
to other bed or to CO2 accumulator. 

= Total flow (lb/hr) through cabin condenser heat 

• CondenIlte f removed (lb/hr) in cabin condenser heat 
exchanger . 


„/«) = Splitting ratio in splitter (ratio of flow to SAWD 

to total secondary flow from cabin). 

The following are input 1 0 * n , d ° ut ^ i^TenTiM^ heat Exchanger): 
tion of component number 5 (mam cabin sensmie 


Input : 
STEADY 
KSYPAS 


VV(2,104) 
VV( 2 ,87) 

Output : 


= Logical variable denoting if steady state calcu- 
1 stinns are now being performed or not. 

= Integer variable which counts the number of passes 
of the computer program through the solution of 

the components. 0 

= Main cabin gas mixture temperature , r 
= Main cabin gas design temperature, f 


w u u v * 

B(1) = Cooling flow to sensible heat exchanger, lbm/hr 

^„„ fo “po"e^ 

Input : 

A(1) = Total primary (gas) flow into humidity control heat 

R(72) . Effective 'speci f i c heat (Cp) of 9« flow.into 

humidity control heat exchanger, Btu/lb-F 
B(1) = Water cooling flow into humidity control 

CPA = Speci fic r ’heat /h of total air f Jow entering the 

CPA humidity control heat exchanger, Btu/lbm- F 
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3.6.1 GP0LY1 (Continued) 



Output : 

r ( 66) = Calculated overall heat transfer coefficient of 

humidity control heat exchanger, Btu/hr-°F 

The following are input to and output from GP0LY1 before the solu- 
tion of component number 10 (water separator for cabin condenser 
heat exchanger). 


Input : 

A ( 7 ) 

A( 1) 
A(5) 

A ( 8 ) 
A(6) 
CPONV 


= Condensable entrained liquid flow, lb/hr 
= Total primary flow into water separator, lb/hr 
= Non-condensables flow into water separator, lb/hr 
= Non-condensables specific heat, Btu/lb-°F 
= Condensable vapor flow into water separator, lb/hr 
= Specific heat of condensable vapor flow into 
separator, Btu/lb-°F 


Output : 

R(67) 

A(l) 

CPA 

R(65) 

A(7) 


= Condensable entrained liquid flow, lb/hr 
= Total flow into water separator after removal of 
entrained liquid, lbm/hr 

= Specific heat of total flow into water separator, 
Btu/lbm-°F 

= Power added to water separator, Btu/hr 
= New entrained liquid flow which is set to 0.0 
lb/hr. 


The following are input to and output from GP0LY1 before the solu- 
tion of component number 14 (splitter for gas flow from main cabin 
to SAWD bed #1 or SAWD bed #2). This section of GP0LY1 is where all 
the control logic for the operation of the SAWD carbon dioxide 
removal system is performed. 


Input : 

STEADY 

TIMEMX 

TIME 
DTI ME 
RHAVG 

TABOLD 

VV( 2,89) 
VV( 21,2) 
VV(21,85) 


= Logical variable denoting if steady state calcula- 
tions are now being performed or not. 

= Length of time the transient case is to run, 

seconds . 

= Length of time into transient, seconds. 

= Time step between interations, seconds. 

= Average relative humidity (decimal form) in cabin 
during absorption cycle. 

= Length of absorption cycle before present itera- 
tion, seconds. 

= Cabin gas relative humidity (decimal form). 

= Temperature of flow exiting SAWD bed #2, °F 
= Exit temp, from SAWD bed to end desorption, F. 
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3.6.1 GP0LY1 (Continued) 


Input : 
VV(2,183) 
KK(21 ,16) 
TBSET 

TEXCH6 

KK( 16 ,16) 
TREND 


VV(16,2) 

A(l) 

RHOA 
VV( 2,20) 

V(IS+I) 


V( IS+6) 
V ( I S +7 ) 


= SAWD absorption cycle time curve multiplication 
= Absorb/desorb flag for SAWD bed #2 (0-absorb, 

= Length'of time during which the bed which jb« other 

finished desorbing bleeds its pressure 

= Maxi mum C length of time (seconds) dur ^ g ^^ e p he 
bed which just finished desorbing will transfer 
energy to the other bed if the other bed does not 

, Absfrb/desorb e flag a for SAWD bed #1 (0=absorb, 

T^^Mtiirp nf SAWD bed receiving energy at which 
‘ S^rans?er W'he other SAWD bed is termin- 

: S'i:r«:'3«s s,K ... - - 

: silrx/iri.r.r v>x«; 

= Thi^statemenlTdenotes « ri “ ! 

lated Ob/hO^nclSde: °*ysen, ‘ nitr09en ’ 

- ;-i" ; whr. 

: Condensable entrained liquid flow from cabin 

. AnonHari/ f 1 OW . lb/hr . 


Output : 

A(I) 

WSTM 

CFM 
A( 1 ) 

VV(14,65) 
V V ( 17 ,65) 
VV ( 22 ,65) 
VV(24,65) 


KX.'K,lsr;; S 

low of steam needed for a particular desorption 

vcle, lb/hr. -a, . 

: low throuqh SAWD bed fan, ft^/min. 

rotal flow 9 to humidity control splitter from cabl , 

b/hr. now 


atio of 
pi i tter 
atio of 
pi i tter 
: at i o of 
iplitter 
tatio of 
splitter 


secondary 

#14. 

secondary 

#17. 

secondary 

# 22 . 

secondary 

#24. 


exit 


to inlet flow for 


exit flow to inlet flow for 
exit flow to inlet flow for 
exit flow to inlet flow for 
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3.6.1 GP0LY1 (Continued) 
VV( 26,65) 

VV(30,65) 


KK( 17 ,8) 
KK(22,8) 
KK(24,8) 
KK( 26 ,8) 
KK( 29 ,8) 


Ratio of secondary exit flow to inlet flow for 

Ratio^of secondary exit flow to inlet flow for 

Numbe^of ^component to be solved after component 

Number of component to be solved after component 

Number of component to be solved after component 

Su*er of component to be solved after component 

Number of component to be solved after component 

#29. 


inlet): 
Input : 
BLEED 
KK( 16 ,16) 


Flag is true if SAWD cycle is in pressure bleed 

Absorb/desorb flag for SAWD bed #1 (0=absorb, 
l=desorb) . 


Output : 

B(I) 


The following are input to and o„tP ; t f™ ' ^ Sj 
"e ’pressure bleed to bed «. 


Input : 

KK(16,16) 

VV( 3 ,23) 

KK( 16 ,19) 

V(IV-KJ) 

R(4) 

R ( 2 ) 

R(5) 

R(6) 

R(9) 


Absorb/desorb flag for SAWD bed #1 (0=absorb, 
Outlet"^ ressure from humidity control splitter on 
Number 3 of SMO bed sfg^ts'^d in calculations in 

SAWD bed #1* . ; n cauq bed #1, ft^ 

oSilerp U r«sSre 5 Tbed #1 for bleeddown, psia 

i Molecular^ weigh^of non-condensables, lb/mol. 
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GP0LY1 (Continued) 


Output : 
PNEW 


= Inlet pressure to SAWD bed #1, psia 


The following are input to and output from t0 

solution of component number 20 (gas mixer steam u 

SAWD bed #2 inlet): 


Input : 

BLEED 

KK(21,16) 


= Flag is true if SAWD cycle is in pressure bleed 

= Absorb/desorb flag for SAWD bed #2 (0=absorb, 
l=desorb). 


Output : 
B ( I ) 


= Secondary side 
zero during bed 


inlet flows which are set equal to 
#2 pressure bleed down to bed #1. 


The following are input to and output ^« # « P “ LV ; e b 4 f ° s re c ‘u U nte<i 
^^ss2 f re C STn dJrl^IrSure bl— to heO ». 


Input : 

KK( 21 ,16) 

KK(21,19) 

V(IV-hJ) 
VV( 3,23) 

R ( 4 ) 

R(2) 

R(5) 

R(6) 

R(9) 


= Absorb/desorb flag for SAWD bed #2 (0 absorb, 

= Number P of * SAWD bed segments used in calculations in 


SAWD bed #2* cm in u a/ h #9 ft 3 

Void volume of segment j in SAWD bed #2. 
Outlet pressure from humidity control splitter 
secondary side (to beds), P sia * 

Outlet pressure at bed #2 for K b l ee ^ d0 ^ ? pSia 
Temperature of fluid exiting bed #2, psia 
Non-condensable flow from bed #2, 1b/hr 
Condensable vapor flow from bed #2, lb/hr 
Molecular weight of non-condensables, lb/mol. 


Output : 

PNEW = Inlet pressure to SAWD bed #2, psia. 

the energy transfer or bleed portions of its cycle. 
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3.6.1 GP0LY1 (Continued) 


Input : 

PREHET = Flag is true if SAWD is in energy transfer portion 

of its cycle. 

BLEED = Flag is true if SAWD is in bleed portion of its 

cycle . 

Output : 

B ( I ) = Flows into secondary side of mixer. 

The following are input to and output from GP0LY1 before the 
solution of component 27 (water supply tank for steam generator). 
Here is calculated the water flow to be used for desorbing a SAWD 
bed . 


Input : 

VV(2,105) 
VV( 32,72) 
VV(30,65) 

STEADY 

R(100) 

KK ( 16,16) 

VV(16,2) 
VV( 16 ,85) 
KK(21 ,16) 

BLEED 

PREHET 

V V ( 2 1 ,2) 
VV (21 ,85) 


= Cabin gas mixture total pressure, psia 
= Pressure in carbon dioxide accumulator, psia 
= Split ratio in component 30 (splitter for flow to 
CO 2 accumulator or to cabin). 

= Logical variable denoting if steady state calcu- 
lations are now being performed or not. 

= Steam flow, pph. 

= Absorb/desorb flag for SAWD bed H (0=absorb, 
l=desorb) 

= Temp exiting SAWD bed #1, °F 
= Exit temperature to end desorption, °F 
= Absorb/desorb flag for SAWD bed #2 (0=absorb, 
l=desorb) . 

= Flag is true if SAWD is in bleed portion of cycle. 
= Flag is true if SAWD is in energy transfer portion 
of cycle. 

= Temp exiting SAWD bed #2, °F 
= Exit temperature to end desorption, °F 


Output : 

R ( 1 ) = Total flow from water tank, lb/hr 

The following are input to and output from GP0LY1 before the 
solution of component 30 (splitter for desorbed gas - to cabin or 
CO 2 accumulator). Here is the control logic to simulate the flow 
sensor sensing flow and to direct the desorbed flow to the accumu- 
lator when the flow reaches the setpoint. 
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3.6.1 GP0LY1 (Continued) 


3.6.2 


Input : 

RHOA 
A( 1) 

KK (16,16) 
KK( 21 ,16) 


Output : 

R(65) 


= Density of flow to splitter, lbm/ft 
= Total flow to splitter, lb/hr 
= Absorb/desorb flag for SAWD bed #1 (0-absor 

l=desorb) . 

= Absorb/desorb flag for SAWD bed 
l=desorb) . 


#2 (0=absorb, 


= Split ratio in splitter which is ratio of flow 
directed to cabin to total entenng flow. 


The following are W /.X - 1 oJK 

ro’^^fneT^now 1..U the - 

the flow is to go overboard or to another system ror 
the carbon dioxide. 


Input : 

STEADY 

R(72) 

KK(32,16) 

A ( 1 ) 

R(71) 

R(70) 


Logical variable denoting if steady state calcu- 
lations are now being performed or not. 

Prpssure in accumulator, psia* , 

Carbon dioxide delivery (0=overboard , 1=C0 2 reduc- 
tion). . , < i 

Total flow to C0 2 accumulator, lb/nr 
C0 2 accumulator tank volume, ft 
C0 2 temperature, °F 


= Total flow out of C0 2 accumulator, lb/hr 


Output : 

R(l) 

GP0LY2 

The GP0LY2 subroutine in ESCM performs "Y*™ 1 . ^waie^ vapo^i^lhe 
the rate of change in oxygen, dar k°" dl “*’ • loadings in P both SAWD 

cabin, to compute the water and card “ . t l water 9 use d for SAWD 

rd S ieso:%)r y anrthVVo\a7?onrsafe%^d7rom the cabin air. 

Input and output f tte “r 1 ^ "a"" 

£ oVt r heTub 0 rL»^eT:errd 0r , k nthe appropriate COMMON b,oc k s. 
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3.6.2 GP0LY2 (Continued) 

The following are input and output from GP0LY2 after component 1 
(crew) is analyzed: 


Input : 

KK(21,16) 

VV( 16 ,72) 
VV( 16 ,73) 
VV(24,65) 


VV(21 ,72) 

VV(21 ,73) 

R(70) 

VV(10,67) 

R(68) 

R(69) 


Output : 

VV(2,137) 

VV(2,175) 

VV(2,177) 


= SAWD bed #2 absorb/desorb flag: 0=absorb, 

l=desorb. 

= Total SAWD bed #1 CO p absorption rate, lbm/hr. 

= Total SAWD bed #1 HgO evaporation rate, lbm/hr. 

= SAWD exit flow control valve 1.0 = flow goes to 
opposite bed for energy transfer; 0.0 = flow goes 
to accumulator. See Figure 4. 

= Total SAWD bed #2 CO 2 absorption rate, lbm/hr. 

= Total SAWD bed #2 water evaporation rate, lbm/hr. 
= Total crew water vapor generation rate, lbm/hr 
= Rate of condensate removal from cabin, lbm/hr 
= Total crew oxygen usage rate, lbm/hr 
= Total crew carbon dioxide generation rate, lbm/hr 


= Net water vapor flow into cabin, lbm/hr 
= Net oxygen addition rate to cabin, lbm/hr 
= Net carbon dioxide addition rate to cabin, lbm/hr 


The following are input and output for GP0LY2 after component 16 
(bed #1) or 21 (bed #2) are analyzed. Here bed loadings are calcu- 
lated, and during steady state, the proper air temperature and mass 
flows of vapor and entrained liquid are computed. 


Input : 


R ( 78) 

R (80) 

R (81 ) 
A(2) 

A ( 5 ) 

A ( 6 ) 

A ( 7 ) 

A(8) 

STEADY 


= Total dry bed #l/#2 resin weight, lbm 
= Total water in bed #l/#2 resin, lbm 
= Total carbon dioxide in bed #l/#2 resin, lbm 
= Bed #1 or bed #2 inlet flow temperature, °F 

= Bed #1 or bed #2 inlet non-condensable flow, lbm/hr 

= Bed #1 or bed #2 inlet water vapor flow, lbm/hr 

= Bed #1 or bed #2 inlet entrained liquid flow, 

lbm/hr 

= Specific heat of non-condensable flow entering bed 
#1 or bed #2, Btu/lbm-°F 

= Logical variable denoting if steady state calcu- 
lations are now being performed or not. 
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3.6.2 GP0LY2 (Continued) 


Output : 

R ( 82 ) 
R(83) 

R(2) 

R(6) 

R(7) 

Note: 


T jn nr bed #2 water loading, lbm/lbm 

I Total bed ft or bed #2 carbon dioxide loading, 

lbm/lbm of dry resin. 

= Bed #l/#2 exit flow temperature, 

= Bed #l/#2 exit water vapor f ow, bm/l h 
= Bed #l/#2 exit entrained liquid flow, / 

Variables R(2). R(6>. and R(7) are output only during steady 
state calculations. 


. , nH nutDut from GP0LY2 after component 
The following are input to and 0lJ tpu 

27 (the water supply tank) is ana y 


Input : 

R(D 

R(98) 

R ( 99 ) 
DTIME 
STEADY 

VV(10,67) 


Total flow exiting water supply tank, lbm/hr 
Previous total water used from tank, lbm 
Previous total water added to tank, lbm 

Unab“ P den e o“"ng S .f steady state calcu- 
lations are now being performed or not. 

Rate of condensate removal from 
(i.e., rate of water added to tank), / 


Output : 

R( 98) 

R(98) 

No error messages are generated by , GP0LY2. Lower level functions 
and subroutines used by GP0LY2 ar 

ISIS u r uc WRI AN f. 


Updated total water used from tank , lbm 
Updated total water added to tank, lbm 


3.7 


SS Subroutine 

e ARSS subroutine generates a schematic ^^J^^o^such as flow, 

vitalization subsystem with perfomance morrn A1$0 

mperature, and heat load printed near eacn comp^^ Jhis sub _ 

i nted are total crew loadings an^cab ^ ^l ther A RSS is called 

utine is called from sburou > through the printoff 

not in ECLST is specified by the use th 9 q vv(2>184) . Figure 

IZTl ’saSe Va sSomltiJ 9“e^ted by E% in tbe ARSS subnet, ne. 
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SPACE STATION AIR REVITALIZATION SUBSYSTEM SAWO II DEMONSTRATION MOOEL. MISSION TIME: 3015. SEC < 0.84 HR) DATE: 01/26/85 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


X 


5 

X 

CABIN 

(2) 

CREW 

(1) 


FROM CABIN 

6 

X 







M= 9357. 

7 

X 

VOLUME 

= 

8000. CU-FT 

NO OF MEN = 

3 


T= 69.5 

8 

X 

AIR MASS 

= 

594.03 LDM 

TOTAL Q = 

1423.5 

B/HR 

A 

4 

X 

TEMP 

= 

69.5 F 

METABOLIC Q = 

474.5 

B/HP/MAN 

A 1 

10 

X 

DEW PT 

= 

51.2 F 

SENSIBLE = 

330.2 

B/HR/MAN 

A 1 

n 

X 

TOTAL PRESS 

= 

14.69 PSIA 

LATENT = 

144.3 

B/HR/MAN 

AAAAAAAAAl 

12 

X 

02 PRESS 


2.89 PSIA 

02 USAGE = 

0.2320 

PPH 

> | 

13 

X 

C02 PRESS 

= 

2.58 MH-HG 

C02 PR 00 = 

0.2733 

PPH 


14 

X 

GAS LEAKAGE 

= 

0.083 PFH 





15 

V 

02 MAKEUP 

= 

0.000 PPH 





16 

X 

N2 MAKEUP 

= 

0.000 PFH 

FROM 




17 

X 

NON ECLSS Q: 



CABIN 




18 

X 

SENSIBLE 

= 

17065. B/HR 

H= 1329. 



19 

X 

LATENT 

= 

255. B/HR 

T= 69.5 


• 


20 

X 

ECLSS Q: 



A 


*T= 56 

.6 

21 

X 

SENSIBLE 

= 

5399. B/Hft 

1 A 

— 



22 

X 

LATENT 

= 

2607. B/HR 

1 A M= 

1388. IREL. HUMD. 

1 


I SENSIBLE 
IHX FAN (4) 


I 

I T = 70.8 


I SENSIBLE 
I HX (5) 


TO CABIN 
M= 9356. 
T = 66.4 

A 

A 

A 


2843. I AAAAAAAAAA l Q=-10105 . | AAAAAAAAAAAAAAAA 


|CFM= 2100. I 


I 


I 


X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 

44 X 

45 X 

46 X 

47 X 

48 X 
X 


REL HUMIOITT = 51.82 PCT 


FROM CABIN 
M= 56.92 
T= 69.5 
A 

A I 
A I 
A V 

A M= 0.0 


V A T= 73.3 I HX (7) 

( 31 ) AAAAAAAAAA I Q=-11026 , 


I T = 50.2 


M= 1344. 


M= 950. 
T= 60.0 


M- 1385. 
T = 50.2 


I REL. HUMD. I 
I HX FAN 112)1 


TO CABIN 
M= 1387. 
T= 51.6 

A 

A A 

A 1 

A I 
A I 


I AAAAAAAAAt 9 )AA AAA AAAAAAAAAA! 11 1AAAAAAAAAAAAAA I Q = 478. I AAA I 


I 


M= 59.17 
T=156.2 


•Ms 950. 
•Ts 45.0 


-> A A 

I A A 

I A I H20 SEP (10)1 A Ms 41.6 

V AA|Q = 13. I AAA T= 51.5 

ICFM = 9. | 


M > H20 INTERFACE: 

UWUWWWUWUUW MS 2.9 


|CFM= 300. 


M ADDED = 1.601 


ISAWD FAN (19)| < 

|Q =171. | AAAAAAAAAA AAA A AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA A 


ICFM= 13. 


M= 1.65- 
T=284 .5 I 


AMINE 


P= 22.77 
— M= 0.0 
I T=110.4 


M= 0.0 


A A 
A I 
A I 


( 14 )AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA( 15 1SSSSSSSSSS I BED 1 (16 ) ICCCCCCCCCCI 17 )AAAAAAAAA< 18 > 


H20 INTERFACE 
M USED = 1.830 

W 
I W 

I W M= 1.65 
V U T= 70.0 


l PUMP 
1128) 


I T= 90.71 STEAM I T = 284.5 


IQ = 1677. I 
|MC02= 0.1951 
|MH20= 2.9241 


TO CABIN 
M= 0.0 
T = 78.5 

A 

A A 
I A 

I A 
I A 
A M= 


| KW*»." JWW.JWW | GEN ( 29 ) I SSSSSSSS( 25/26 )SSSSSSSSSSSSSSSSSSSSSSSS5SSSSS( 23/24 ICCCCCCCC )A( CCCC( 30 JCCCCCCCCl 


I C02 I 

IACCUM (32)1 


A |Q= 34.1 


I Q= 1860.1 


I 

I M= 0.379 I 


50 
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Figure 6 

Sample Schematic Generated By ESCM 
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ARSS Subroutine (Continued) 


The following are input into S“ bro h u ^ n ® *ept ^hos^ wU^asterisks . 
no manipulations performed ojJi hem . they are output. The 

are ’ n P ut - 

not. outDUt. 


I DATE 
TIME 

The following 


= Date of transient run 
= Time into transient. 


(input and output), 
seconds (input and output). 


inputs/outputs relate to the crew (component 1): 


VV(1,65) 
V( IV+2) 

V ( IV+3) 

V ( I V +4 ) 

V ( I V +5 ) 
KK( 1 ,16) 


= Total metabolic heat from crew, Btu/hr 
= Sensible heat per crewman, Btu/hr 
= Latent heat per crewman, Btu/hr 
= Total oxygen use rate, Ib/hr 
= Total CO 2 generation rate, lb/nr 
= Number of crewmen 


The following inputs/outputs relate to the main cabin 


VV(2,1) 

V ( I P +2 ) 
V( IV+2) 
*V(IV+25) 


V( IV+30) 

V ( I V +3 4 ) 

V( IV +36) 

V (IV +39 ) 
V(IV+40) 
V( IV+42) 

V ( I V +58 ) 
**V( IV +64) 
**V( IV+72) 

V( IVT+1) 
V( IVEX+6) 
V( IVEX+7) 


Total flow out of cabin to sensible heat exchanger 
Temperature of flow exiting cabin and flowing to 

Equipment^ lighting and miscellaneous heat load 

added to cabin gas, Btu/hr 
Cabin qas relative humidity, decimal 
fornn This is multipled by 100 to put it 
in percent form before it is output. 

Cabin gas oxygen presure, psia 
dew point, °F 

carbon dioxide pressure, mmHg 
mixture total mass, lbs 
mixture temperature, F 
uuu. yao mixture total presure, psia 
Outboard leakage rate, 1b / h ^ 

Water vapor addition rate, lb/nr 

_ i_ i ^ £1 A 1 .I 4 nf rt r AT 


Cabin 
Cabi n 
Cabi n 
Cabi n 
Cabin 


gas 

gas 

gas 

gas 

gas 




lb/hr 3 

= Cabin free volume, ft 
= Diluent addition rate, lb/hr 
= Carbon dioxide addition rate. 


lb/hr 
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ARSS Subroutine (Continued) 

The following inputs/outputs relate to the splitter (component 3) 
for flow from the cabin: 


VV (3,1) = Total flow from splitter (3) to gas mixer (31) then 

to humidity control heat exchanger, lb/hr 
V ( I P +2 ) = Fluid temperature exiting splitter primary side 

(3), °F 

V ( IS+1 ) = Total flow from splitter (3) to SAWD beds, lb/hr 

V ( I S +2 ) = Fluid temperature exiting splitter (3) to SAWD 

beds, °F 

The following inputs/outputs relate to the main cabin sensible heat 
exchanger fan (component 4): 


VV( 4,2) = Temperature of fluid exiting main cabin sensible 

heat exchanger fan, °F 

V(IV+1) = Heat added to gas stream at fan, Btu/hr 

V( IV+12) = Cabin fan volumetric flow rate, cfm 


The following inputs/outputs relate to the sensible heat exchanger 
(component 5): 


VV (5,1) 

V ( I P +2 ) 

V ( I S +2 ) 
V(IV+1) 
VV (6,1) 
W (6,2) 


= Total flow from sensible heat exchanger (5) (prim- 
ary side - gas side), lb/hr 
= Temperature of gas exiting sensible heat exchanger 
(primary side), °F 

= Temperature of fluid exiting sensible heat ex- 
changer (secondary side), °F 
= Heat transferred across heat exchanger between 
secondary and primary sides, Btu/hr 
= Coolant flow to cabin sensible heat exchanger (5), 
lb/hr 

= Temperature of coolant to cabin sensible heat 
exchanger, °F 


The following inputs/outputs relate to the main cabin condensing 
heat exchanger (component 7, also known as humidity control heat 
exchanger) : 


VV (7,2) 

V ( I S +2 ) 

V(IV+1) 
W(8,l) 
VV( 8,2) 


= Temperature (°F) of fluid exiting main cabin 
condensing heat exchanger 7 (humidity control heat 
exchanger), °F 

= Temperature (°F) of fluid exiting main cabin 
condensing heat exchanger 7, secondary or coolant 
side . 

= Heat transferred across condensing heat exchanger 
7, Btu/hr 

= Coolant flow to main cabin condensing heat ex- 
changer 7, lb/hr 

= Temperature of coolant to cabin condensing heat 
exchanger, °F 
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ARSS Subroutine (Continued) 


VV (9,1) 


Total flow exiting splitter (9) and flowing to gas 
mixer (11)» lb/hr 


+ c /rt.itnuts relate to the water separator (com- 
The following inputs/outputs relate 

ponent 10): 


VV(10,1) 

V( IP+2) 

V( IV+1) 
V ( I V +3 ) 


, Total flow from water separator (10) to gas mixer 

. Temperature of fluid exiting water separator ,10, 
and going to gas mixer (ID. 

= Condensate * removed in water separator, lb/hr 


/ * 4 -c roiatP to the gas mixer (component 11) 

for fio^t^condense^heat' ex changer fan: 


VV(ll.l) 
VV( IP+2) 


. Total flow from gas mixer (11) to condenser heat 

. Temperature^ ' fluid^ flowi ng fr» gas mixer (11) to 
condenser heat exchanger fan, r 


tc/m.tnuts relate to the condenser heat exchanger 
The following inputs/outputs relate 

fan: 


VV( 12 ,1) 

V( IP+2) 

V(IV+1) 

V ( I V +1 2 ) 


= Total flow exiting condenser heat exchanger fan 

. Temperature of fluid flowing from condenser heat 

- Total "heat^added^ to gas stream Btu/hr 
= Volumetric flow through fan (12), ctm 


The following inputs/outputs relate to the splitter (component 14, 
for flow to SAWD beds: 


VV( 14,1) 

V(IS+1) 

V(IV+1) 


= Total flow to SAWD bed #1 from splitter (M), lb/ hr 

= Total flow to SAWt J r # M 4 )°f 0 r P flow e to SAWD beds. 
= Split ratio in splitter (14) tor now 


* / +m.tc rplate to the gas mixer (component 15) 
The following inputs/outputs rela 

for flow to SAWD bed #1: 


VV(15 ,1) 
V ( I P +2 ) 


Total flow to SAWD bed *1 from gas mixer (15), 

Temperature of fluid exiting gas mixer (15) and 
going to SAWD bed #1, F 
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ARSS Subroutine (Continued) 

The following inputs/outputs relate to SAWD bed #1 (component 16): 


VV(16,1) 

V ( I P +2 ) 

V ( IP+4) 

V ( I V +6 ) 

**KK(16 ,16) 
**V( IV+9) 

V ( I V +16) 

V( IV+17) 


= Total flow exiting SAWD bed #1, lb/hr 
= Temperature of fluid exiting SAWD bed #1, °F 
= Outlet pressure at SAWD bed #1, psia 
= Gas stream total sensible heat change in bed #1, 
Btu/hr 

= Absorb/desorb flag (0=absorb, l=desorb). 

= Total water evaporation/condensation rate (pph) 

= Total water in SAWD bed #1 at start of transient, 

^s . 

= Total C0 2 in SAWD bed #1 at start of transient, lbs 


The following inputs/outputs relate to the splitter (component 17) 
for flow from SAWD bed #1: 


VV (17,1) = Total flow exiting splitter (17) which goes to the 

condenser heat exchanger, lb/hr 

V(is+1) = Total flow exiting splitter (17) which goes to 

C0 2 recovery, Ib/hr 

V(IV+1) = Split ratio in splitter (17). 

The following inputs/outputs relate to the SAWD bed fan (component 
19): 

VV( 19 , 1) = Total flow through SAWD bed fan (19), lbm/hr 

V ( I P +2 ) = Temperature of flow exiting SAWD bed fan, °F 

V(iv+1) = Total heat added to gas stream, Btu/hr 

V ( I V +1 2 ) = SAWD bed fan volumetric flow rate, cfm 

The following inputs/outputs relate to the gas mixer (component 20) 
for flow to SAWD bed # 2 : 


VV(20,1) = Total flow from gas mixer (20) which sends steam or 

cabin gas to SAWD bed #2 inlet, lb/hr ^ 

V(IV+2) = Temperature of flow from gas mixer (20), °F 

The following inputs/outputs relate to SAWD bed #2 (component 21): 


VV(21,1) 

V ( I P +2 ) 

V( IP+4) 

V ( IV+6) 
**KK( 21 ,16) 
**V( IV+9) 

V( IV+16) 

V( IV+17) 


= Total flow exiting SAWD bed #2, lb/hr 
= Temperature of fluid exiting SAWD bed #2, °F 
= SAWD bed #2 outlet pressure, psia 
= Gas stream total sensible heat change, Btu/hr 
= Absorb/desorb flag (0=absorb, l=desorb). 

= Total water evaporation/condensation rate, lb/hr 
= Total water in SAWD bed #2 at start of transient, 
lbs 

= Total C0 2 in SAWD bed #2 at start of transient, lbs 
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ARSS Subroutine (Continued) 

The following inputs/outputs relate to the splitter (component 22) 
for flow from SAWD bed #2: 

VV(22,1) = Total flow exiting splitter (22) which goes to 

condenser heat exchanger, lb/hr 

V ( I S +1 ) = Total flow exiting splitter (22) which goes to 

CO 2 recovery, lbm/hr 

V(IV+1) = Split ratio in splitter (22) 

The following inputs/outputs relate to the splitter (component 24) 
for flow from gas mixer (23): 


VV(24,1) = Total flow exiting splitter (24) to CO 2 accumu- 
lator, lb/hr . 

V ( 1 s +1 ) = Total flow exiting splitter (24) to gas mixer (2b), 

lb/hr 

The following inputs/outputs relate to the splitter (component 26) 
for flow from gas mixer (25): 


VV( 26 ,1 ) = Tota flow of steam from splitter (26) to SAWD bed 

V(IS+1) = Total flow of steam from splitter (26) to SAWD bed 

#2, lb/hr 


The following inputs/outputs relate to the water supply tank (com- 
ponent 27): 


VV( 27 ,1 ) 

V(IP+2) 
V( IEX+1 ) 
V ( I E X +2 ) 


= Total flow from water supply tank (27) to steam 
generator, lb/hr 

= Temperature of flow from water supply tank, F 
= Total water used, lb/hr 
= Total water added, lb/hr 


The following inputs/outputs relate to the water supply pump (com- 
ponent 28): 


VV( 28,2) 
V(IV+1) 


= Temperature of water exiting water supply pump (28) 
for steam generator, °F 

= Total heat added to water in pump (28), Btu/hr 


The following inputs/outputs relate to the steam generator (com- 
ponent 29): 


VV(29 2) = Temperature of steam exiting steam generator 

(29) °F 

*V(IV+4) = Electric energy required by steam generator, watts. 

It is multipled by 3.413 to convert it to Btu/hr 
before it becomes output. 
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ARSS Subroutine (Continued) 

The following inputs/outputs relate to the splitter 
for flow from gas mixer (24): 


(component 30) 


VV(30,1) 

V(IS+1) 

V(IS+2) 

V(IV+1) 


= Total flow exiting splitter (30) to C0 2 accumu- 

= Total "flow exiting splitter (30) to cabin, lb/hr 
= Temperature of fluid exiting splitter (30) to 

= Ratio*of flow directed to cabin to total inlet flow 
in splitter (30). 


The following inputs/outpots relate to the gas nrixer (component 31) 
for flow to condensing heat exchanger (/)• 


VV(31,1) 
V ( I P +2 ) 


= Total flow exiting gas mixer (31) to humidity 
control heat exchanger (7), lb/hr 
= Temperature of fluid exiting gas mixer (31) to 
humidity control heat exchanger, F 


The following inputs/outputs relate to the C0 2 accumulator (com- 
ponent 32): 


VV(32,1) 

V ( I P +2 ) 
V( IV+5) 
V( IV+8) 


= Total flow exiting C0 2 accumulator (32) to over- 
board or C0 2 reduction, lb/hr 
= Temperature of fluid exiting C0 2 accumu ator, F 
= Total fluid weight in C0 2 accumulator, lbm 
= C0 2 accumulator pressure, psia. 


The following input relates to gas mixer (component 18). 


**VV( 18 ,6) 


= Condensable vapor flow exiting gas mixer (18), 
lb/hr 


The following input relates to the splitter (component 3) for flow 
from the cabin: 


**VV(3 ,25) 


= Condensable vapor flow exiting splitter 
(3) and flowing to SAWD system. 


The following are output only of ARSS: 


TIMEM 


= Time into transient in minutes. It is calculated 
by dividing the input TIME (seconds) by 60. 


TIMEH 

QMET 


ime into transient in hours. It is calculated by 
lividing the input TIME (seconds) by 3600. 

‘otal metabolic heat per crewman, Btu/hr. it is 

_ . . _ • L 1 - 1 if a 


1 rt-f 1- ho 


1 . j ^ + r 


per crewman* 
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ARSS Subroutine (Continued) 


3.8 


= Latent heat added to cabin by non-ECLS and non- 
metabolic sources, ie, vehicle loads, Btu/nr. 
= 9.0 = water separator power, Btu/hr 
= Latent heat released due to condensation during 
desorption or sensible heat used for evaporation 
during absorption in SAWD bed #1, Btu/hr 
= Latent heat released due to condensation during 
desorption or sensible heat used for evaporation 
during absorption in SAWD bed #2, Btu/hr 
= Total power required by ECLS components, Btu/hr 
= Total latent heat added to cabin from ECLS sources, 
i.e., the SAWD system, Btu/hr 

The following are outputs which cause A's, C's or s s . w ^ch rep- 
resent air, carbon dioxide, and steam flows, respective y 
printed at the proper places and times on the output schematic. 


QNELAT 

V1076 

Q16 


Q21 


QESEN 

QELAT 


LI 

L2 

L3 


L4 

L5 

L6 


L7 
L8 

IR45 Subroutines 


Flow stream into SAWD bed #1 

Flow stream out of SAWD bed #1 to splitter (17) 

Flow stream between splitter (17) and gas mixer 
(23) which ultimately goes to bed #2, the cabin, 
or the C0 2 accumulator. 

Flow stream into SAWD bed #2 

Flow stream out of SAWD bed #2 to splitter (22) 

Flow stream between splitter (22) and gas mixer 
(23) which ultimately goes to bed #1, the cabin, 
or the C0 2 accumulator. 

Flow stream between splitter (24) and the C0 2 
accumulator. 

: Flow stream between splitter (30) and the cabin. 


the transient performance 


of a 
exit 
be 


This qroup of subroutines computes 

SAWD canister which includes the inlet head ® r ’ a 

header, and an exit temperature sensor. The bed itself 
divided into as many as five segments. 

The computational sequence for subroutine IR45 is shown in 1 Figure 7. 
At the start of the subroutine, values of resin weight, carbon 
dioxide weight, water weight, canister weight, ™id volura, and fo ™ 
weight are computed for each segment of the bed • This is done only 
the first time that the IR45 subroutine is called in ISC*. If 
steady state analyses are to be done, subroutine STAEDY is called, 
if not, the program begins the transient IR45 analysis with a call 
to subroutine HEADER. In that subroutine, the mixed temperature and 
composition of the gases in the inlet header computed F om 

there, the gases enter the first segment of the bed. Subroutine 
BALNCE performs a mass, energy, and pressure balance on the gas 


and 
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Figure 7 

Subroutine IR45 Flow Diagram 
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IR45 Subroutines 


SfwcVcaVs’X^ 

SSSs r 1 5HS « 

'riS o? V" til n a rsj n VegJnt a^mpuSd *?S 

SSis ^avl SSB 2 

the temperature and co mp°sition of the ga y , va i ues 

The last subroutine called is UPDAYT. mere, cne / f th 

sS temperature „ compete,. 
IR45 is the main subroutine which simulates a ^AWD bed, the sub^ 

sa: ^ and 

output to this subroutine are pesented in the following. 

A array data and R array data up to R<64) [R<1> to R(64)] are always 
the same Only the component to which they refer changes. 

A Array Data (Input/Output) 

The A array [A( 1) to A( 19)] contains the primary source flow data 
(upstreart olmponent flo, [data which enters the prinmry .,*•_* * 

rrpSf?c’?nVtrf e or^ray'UtioV^ be found in the 0189 
Manual (Table 2-1, page 2-12). 


R Array Data (Input/Output) 


R array locations R(l) to R(19) contain the exit f low c f °r the 

involves exit rather than source flaws. The specific inputs for 
r ( 1 ) to R( 19) may be found in the G189 Manual (page 2-12 

D/?ni tn r f 62 1 are not used in the IR45 subroutine. R(63) and R(64) 
In the IR45 subroutine they refer to SAWD bed #1 or #2. 


R(63) 

R(64) 


= Exterior surface temperature (°F) 

= Conductance between source (bed) and exte 
(Btu/hr-°F ) 


ri or 
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IR45 Subroutines (Continued) 

R(65) to R(95) 

These are subroutine dependent data. Inputs to these positions 
depend upon the subroutine being used (in this case IR45). In IR45 
the R(65) to R ( 95 ) refer to the following in the SAWD beds. 


R(65) 
R(66) 
R(67) 
R(68) 
R ( 69) 
R(70) 
R ( 71 ) 
R ( 72 ) 
R(73) 
R(74) 


= Heat of CO 2 reaction with resin bed (Btu/lb) 

= Specific heat of dry resin bed (Btu/lb- F) 

= Density of dry resin bed (lb/ft ) 

= Resin bed void volume fraction (ft /ft ) 

= Time required to complete steam desorb (sec) 

= Gas stream total sensible heat change (Btu/hr) 

= Gas stream total latent heat change (Btu/hr) 

= Carbon dioxide absorption rate (lb/hr) 

= Water absorption rate (lb/hr) 

= Basepoint absorption/desorption rate (pph of 

C02/1 b of dry resin/psi of CO 2 partial pressure 

difference). , 

= Average steady state desorption steam rate (Ib/nrj 
= Steady state carbon dioxide removal efficiency. 

= Average steady state bed temperature (°F) 

= Total dry resin bed weight (lbm) 

= Total canister weight (lbm) 

= Total water in bed (lbm) 

= Total carbon dioxide in bed (lbm) 

= Total bed water loading (lbm H20/1 b bed) 

= Total bed carbon dioxide loading (lb CO 2 /I o bed) 

= Specific heat of canister (Btu/lbm- 0 F) 

= Temperature to end desorption (°F) 

= Time at which absorption is to end (sec) 

= Outlet temperature time lag (sec) 

= Bed inlet dew point (°F) 

= Bed exit dew point (°F) . 

= Bed inlet carbon dioxide partial pressure (mmHg) 

= Bed exit carbon dioxide partial pressure (mmHg) 

= Total bed exit flow (ft 3 /min) 

= Total foam in bed (lbm) 

= Specific heat of foam (Btu/lbm- F) 

= are spare positions in the array. 


R(75) 

R ( 76 ) 

R ( 77 ) 

R(78) 

R ( 79) 

R(80) 

R(81) 

R ( 82 ) 

R(83) 

R(84) 

R(85) 

R(86) 

R(87) 

R(88) 

R(89) 

R(90) 

R(91) 

R(92) 

R(93) 

R(94) 

R(95) to R ( 99) 

The following R array positions are used for each segment t that the 
SAWD beds are divided into for calculation purposes: In the follow- 

ing J=99+30*(k-l) and k=l thru number of bed segments used. 


R ( J +1 ) = Same as A(I) array 

R(J+20) = Segment k bed temperature ( F) 

R ( J +2 1 ) = Segment k bed carbon dioxide weight (lbm) 

R(J+22) = Segment k bed water weight (lbm) 

R(J+23) = Segment k bed weight (lbm) 
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4.1.1 


I R45 Subroutines (Continued) . 


R(J+24) 

R( J+25) 

R( J+26) 

R ( J +2 7 ) 
R(J+28) 
R(J+29) and 
R( J+30) 


= Segment k canister weight (lbm) 

= Segment k void volume (ft 3 ) 

= Segment k carbon dioxide loading (lb C02/lb bed) 
= Segment k water loading (lb ^O/lb bed) 

= Segment k foam weight (lbm) 

= Spares. 


PROGRAM USE 


Instructions to access ESCM and to use the program are presented in 
this section. ESCM may be run interactively at a terminal or run by 
submitting batch jobs. Of course the job control language will vary 
from computer system to computer system. Also, several subroutines 
will be different from system to system. The job control language 
and subroutines described in section 4.1 apply to the IBM 3080 
series computer with the MVS 3.8 Operating System. Section 4.2 
describes how the ESCM program must be modified to operate on the 
LaRC PRIME computers under FORTRAN 77. 

Operation Using Hamilton Standard's IBM 3080 Series Computer 
To run ESCM, use the following sequence of steps: 

(1) Access the desired command list: 

(a) ESCMCN CNTL (to run batch job) 

(b) ESCMCL CLIST (to run interactive) 

(2) Access input data set: ESCM2B.DATA 

(3) Make desired changes to ESCM2B.DATA 

(4) For batch jobs, SUBMIT ESCMCN; for interactive jobs, EX ESCMCL. 
Accessing Desired Job Control List 

After logging onto the computer time sharing system, either ESCMCN 

or ESCMCL must be accessed from the library and brought into the 

User's file. The following are the commands to do this: 

LC ESCMCN CNTL S0(G15UARL) RON 
LGCL ESCMCL G15UARL RON 

A listing of each is given in Appendix A. Input output files are 
allocated by these control lists; Table 7 presents these assign- 
ments. The program has been coded in such a fashion as to facili- 
tate changes in 1/0 assignments to conform to individual computer 
installations. This was accomplished by using variable unit names 
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Table 7 


FORTRAN 
Vari able 

Current 

FORTRAN 

Logical 

ESCM 

I/O File 

Physical 

Assignments For IBM Computer 

Unit 

Unit 

File 

Unit 


Name 

Number 

Type 

Assignment 


NDA 

1 

Scratch 

These 

Contains card images of all input data cards. Editing and 
processing take place from this file. 

NOB 

2 

Scratch 

four 

files 

may 

use 

Used during Basic Case in conjunction with NDC and NDB to 
process cards. At all other times, contains edited Basic 
Case K/Y array. 

NOC 

3 

Scratch 

Tape, 

Disk 

or 

Drum 

Units 

Used during Basic Case input operations with NDB to process 
K and V Component Data card images. Used to store processed 
S0-4060 plot data at conclusion of 6189A case. 

NDD 

4 

Scratch 

Used during Basic Case with NDB to process Table Data card 
images. May be used for SD-4060 operations. 

NDE 

16 

Output 

Tape 

Output file for G189A plot data to be processed by post 
Processor plot package (Section 10 of Reference 3). 

NTIN 

10 

Input 

Tape 

Input file for the data tape. 

NT0UT 

11 

Output 

Tape 

Output file for the data tape. 

NSTUFF 

31 

Output 

Card 

Punch 

Used to punch card images of the edited simulation data. 



20 

Output 

Terminal 

Output to terminal screen. 
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Accessing Desired Job_Controi^ (Continued) 

URTTEINDA 30)K. The variable 

in all FORTRAN I/O ^f^lTto^he appropriate subprograms in the 
unit names are made available to i ,/linoc the variable unit 

"oln block /FILES/. The ? “IbLOC^ for more infor- 
names are set in the block data subrou^ ^ j$ not set up , , the 

S;rV^on e Tme 3 20 which ?. only used in subroutine ARSS. 


4.1.2 Input Data 


Inout to the program is made through the data set ESCM2B. 
{r 5 da?a sit from Librarian, enter: 


To access 


LG ESCM2B DATA S0(G15UARL) 


L.U LJVyl ~ ' 

A copy of ESCM2B is shown in Table >. B. The ™P u * e * nd up V for the ESCM 
ied in Table 8 are the re conrnjBnded v thgt require d for G189A; 

Reference 3 for more detailed information. 

The input is divided into four 9 ^?Qt S data Card The following sections 
component data, table data and each of these sections 

briefly describe th F V M P the Jser 5gain should refer, to the G189A 
specific for the ESCM, * rom oleteand general description of the 
Manual, Reference 3, for acomple J i pr o V ided to show the general 
innut For convenience. Figure o • * w 
input setup and the case data required. 


4. 1.2.1 Case Data Grou£ 


case Ligt- o w.' w/“r_ 

k a*** AC necessary to specify the 

The case data group cont ^ S me ^ Ch d es^ibe the model in general , and 
overall data storage arrangement, des ^ gpoup als0 contains fluid 
control the execution o blicat i on 0 f Reference 3, an additiona 

input 'has-been adied to th?s group, and it is the first recor 

the group. 


Record 1 


FORMAT cols 


1) Type of printout desired: 


14 


1-4 
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OF 


Table 8 
Input to ESCM 


***» TSO FOREGROUND HARDCOPY ****14:55:21 85030 
DSNAME =TSCG156 . ESCM2B .DATA 

0 o PRINT INPUT DATA AND SCHEMATICS .. . .. 

**** NAMELIST, BASIC CASE 

TAPE 

BASIC 1 50 13 2 YEA NAY 

CA SE ESCM: BASELINE, VOL=8000 FT-3, FC02=2.5, ,275 PPH AVG 

«CASE1 

KPTINVt 1 ) = 1, KPTINVC 2 ) = 5, KPTINV(3)= 10, KPTINV( 4 1=30 , 

KCHOUT= 0, KPUNCH= 0, KPRNT= 0, KRUN= 1, KSTEDY= 0, — 

MAXSLP= 4, MAXSSI= 15. MINSSI= 10, NOLIM=0, 

DTIME=15., START= 0.0, TIMEMX= 66400., 

TMAX= 325.0, TMIN= 0.0, WTMAX= 1.6E6 SEND 


CASE 

GROUP 

DATA 


$ PROP I 

CP<1) = 1.» RH0(1)=62.4, VISC( 1 ) = 2 .42 » WTM( 1 ) = 18. 016 , XK(1)=0.34> 

CP( 2 ) = 0 . 2 , PHO(2)=0. 00045, VISCC 2 ) = 0 . 378, WTM( 2 )=44. 01 , XK( 2 )=0 . 0093, 

CPCONL=l . , CPCONV=0.445,CPCO2=0 .2, CPDIL=0.25, CPOXY=0.22 » CPTC=0 .2, 

GAMGA5=lT4^VISGAS= . 44 »WTMCCN=18 . 016 , WTM0IL=23 . 013»WTMTC=20 . > 

XKGAS=0 . 146 SEND 


ID** 

1 

ARS OEMI 

ID** 

1 

10 FLUID TY! 

ID** 

1 

11 1 

ID** 

1 

12 2 

ID** 

1 

90 

ID** 

1 

100 CREW 

KBAS 

1 

2 1 

NSTR 

1 

0 

VARY 

1 

2 50.0 

VARY 

1 

3 14.7 

VARY 

1 

6 87.8 

VARY 

1 

10 2046. 9 

VARY 

1 

11 7200.7 

VARY 

1 

12 46.7 

KARY 

1 

16 

VARY 

1 

65 

VARY 

1 

66 291.7 

VARY 

1 

'67 175.2 

VARY 

1 

66 

VARY 

1 

69 

VARY 

1 

70 

VARY ‘ 

1 

75" 

VARY 

1 

76 

VARY 

1 

77 

VARY 

1 

78 


T 


«> 70 F > 3 MM HG 


ID** 
KB AS 
MSTR 
VARY 
VARY 
VARY 
VARY 
VARY 
VARY 
VARY 


2 

2 

2 

2 

2 

2 

2 

Z 

2 

2 


3 

6 

10 

11 

12 
21 


MAIN CABIN 
1 25 

151011000 
2 70. 

14.7 
69.7 
2046.6 
7199.4 
46.3 
70. 


5 2 

1 NO TESTS, SS CALCS 

TEMP (F) _ 

► PRESS ( PSIA ) 

H20 VAPOR FLOW (LB/HR) 

OXYGEN FLOW (LB/HR) 

NITROGEN FLOW ( LB/HR ^ _ _ _ , 

C02 FLOW (LB/HR) 

NUMBER OF CREWMEN 

TOTAL METABOLIC HEAT - ALL CREWMEN (BTU/HR 
SENSIBLE HEAT PER CREWMAN ( BTU/HR ) 

""LATENT HEAT PER CREWMAN (BTU/HR) 

TOTAL OXYGEN USE RATE (LB/HR) 

TOTAL C02 GENERATION RATE (LB/HR) 

TOTAL H20 VAFOR GENERATION RATE ( LB/HR ) 
INLET GAS RELATIVE HUMIDITY (DECIMAL) 

INLET GAS DEW POINT (F) 

OUTLET GAS RELATIVE HUMIDITY( DECIMAL ) 
OUTLET GAS DEW PONT (F) 

TOTAL METABOLIC RATE PER CREWMAN (BTU/HR ) 

12 -13 2 

PRI P, 5 RESETS, SPECIFY NET FLOWS 
CABIN GAS MIXTURE: INITIAL TEMP (F) 

4 TOTAL PRESSURE (PSIA) 

H20 VAPOR FLOW (LB/HR) 

OXYGEN FLOW (LB/HR) 

NITROGEN FLCN r LB /HR ) 

C02 FLOW (LB/HR) 

CA3IN GAS MIXTURE: INITIAL TEMP (F) 


14 


START OF 
OF K AND V 
COMPONENT 
DATA 
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Table 8 (Continued) 
Input to ESCM 


VARY 

2 

22 14.7 23 

VARY 

2 

25 13.00 1 

VARY 

2 

29 303.76 i 

VARY 

2 

30 1068.67 

VARY 

2 

31 9.45 J 

VARY 

2 

65 

VARY 

2 

66 17065. 

VARY 

2 

71 

VARY 

2 

72 

VARY 

2 

73 

VARY 

2 

74 

VARY 

2 

87 70.0 

VARY 

2 

88 10.0 

VARY 

2 

89 

VARY 

2 

90 14.7 

VARY 

2 

91 0.5 

VARY 

2 

92 2.9 

VARY 

2 

93 0.5 

VARY 

2 

94 

VARY 

2 

95 

VARY 

2 

96 50.0 

VARY 

2 

97 30.0 

VARY 

2 

98 

VARY 

2 

99 20. 

VARY 

2 

100 

VARY 

2 

101 250. ~ 

VARY 

2 

102 

VARY 

2 

103 

VARY 

2 

104 

VARY ~ 

2 

105 “ 

VARY 

2 

107 

VARY 

2 

108 

VARY 

2 

109 

“ VARY ~ 

"2 

110 ’ 

VARY 

2 

111 

VARY 

2 

112 

VARY 

2 

113 

VARY 

2 

114 

VARY 

2 

115 

VARY 

2 

122 0.0833 

VARY 

2 

123 

VARY' 

2 

124 

VARY 

2 

125 

VARY 

2 

126 

VARY 

2 

127 

VARY 

"2 

128 0.233 ‘ 

VARY 

2 

129 0.0 

VARY 

2 

130 

VARY 

2 

132 70.0 

VARY “ 

~ 2 

135 

VARY 

2 

136 


(PSIA) 

H20 VAPOR FLOW (LB/HR) 

OXYGEN FLOW (LB/HR) 

NITROGEN FLOW (LB/HR) 

CO 2 FLOW { LB/HR) 

TOTAL HEAT ADDED TO CABIN GAS MIXTURE (BTU/HR) 

CABIN HEAT LOAD (BTU/HR) = 5 KW 

HEAT LOSS DUE TO OUTBOARD LEAKAGE (BTU/HR) 

HEAT GAIN DUE TO MASS ADDITIONS (BTU/HR) 

HEAT REQ'D TO FLASH ENTRAINED H20 (BTU/HR) 

FLASH EVAPORATION RATE OF ENTRAINED H20 ( LB/HR ) 
CABIN GAS DESIGN TEMP (F) (4 MEN) 

CABIN GAS OESIGN TEMP TOL (F) (4 MEN) ... 

CABIN GAS RELATIVE HUMIDITY (DECIMAL) 

DESIGN TOTAL PRESSURE (PSIA) 

DESIGN TOTAL PRESSURE TOL (PSIA) 

DESIGN OXYGEN PRESSURE (PSIA) 


CONTINUATION 
OF K AND V 
COMPONENT 
DATA 


VARY 

VARY 

VARY 

VARY 

VARY 

VARY 

VARY 

VARY 

VARY 


137 

139 

17D 

171 

172 

173 

174 

175 

176 


8000. 

70.0 

70.0 

70.0 

70.0 

70.0 


DESIGN OXYGEN FRESSURE TOL (PSIA) 

CABIN GAS OXYGEN PRESSURE (PSIA) 

CABIN GAS NITROGEN PRESSURE (PSIA) 

DESIGN DEW POINT (F) 

DESIGN DEW POINT TOLTF) 

CABIN GAS DEW POINT (F) 

MAX ALLOWABLE CO 2 PRESSURE (MM KG) 

CABIN GAS C02 PRESSURE (MM HG) 

MAX ALLOWABLE TRACE CONTAMINANT LEVEL (PFM) 

CABIN GAS TRACE CONTAMINANT LEVEL (PPM) 

CABIN GAS MIXTURE, TOTAL MASS (LB) 

CABIN GAS MIXTURE, TEMPERATURE (F) _ 

CABIN GAS MIXTURE, TOTAL PRESSURE (PSIA) 

CABIN GAS MIXTURE, NCN-CCNDESABLES MASS (LB) 

CABIN GAS MIXTURE, H20 VAPOR MASS (LB) 

CABIN GAS MIXTURE, ENTRAINED H2Q MASS (LB) 

"CABIN GAS MIXTURE, NON-CONO. SPECIFIC HEAT (BTU/LB-F 
CABIN GAS MIXTURE, NON-CONO. MOL WEIGHT (LB/MOL) 
CABIN GAS MIXTURE, OXYGEN MASS (LB) 

CABIN GAS MIXTURE, NITROGEN MASS (LB) _ _ 

CABIN GAS MIXTURE, C02 MASS ( LB ) 

CABIN GAS MIXTURE, TRACE CONTAMINANTS MASS (LB) 
OUTBOARD LEAKAGE RATE (LB/HR) 

NON-CONDENSABLES LEAKAGE RATE (LB/HR) 

~H20 VAPOR LEAKAGE RATE (LB/HR) 

ENTRAINED H20 LEAKAGE RATE (LB/HR) 

TOTAL MASS ADDITION RATE (LB/HR) 

NON-CONO . ADDITION RATE (LB/HR) __ 

H20 VAPOR ADDITION RATE (LB/HP) 

ENTRAINED H20 ADDITION RATE (LB/HR) 

AVERAGE TEMPERATURE OF MASS ADDITIONS (F) 

ADD COMD VAPOR TEMP 

NET TOTAL FLOW INTO CABIN T LB/HR) 

NET NON-CONO. FLOW INTO CABIN (LB/HR) 

NET H20 VAPOR FLOW INTO CABIN - CALC IN GFOLY2 
CABIN FREE VOL ( FT3 ) 

OXYGEN ADDITION TEMP CF) “ 

NITROGEN ADDITION TEMP (F) 

C02 ADDITION TEMP (F) 

TRACE CONTAMINANTS ADDITION TEMP (F) 

SPECIAL FLOW NO. 1 ADOTICN TEMP (F) 

NET 02 ADDITION (LB/HR) - GPOLY2 CALC 
NET N2 ADDITION (LB/HR) - GPOLY2 CALC 
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Table 8 (Continued) 
Input to ESCM 


VARY 
VARY 
VARY 
VARY 
VARY 
VARY 
VARY 
VARY 
ID** _ 
' KBAS 
NSTR 
VARY 
ID** 
KBAS 
NSTR 
VARY 
VARY 
VARY 
ID** 
KBAS 
NSTR 
' VARY 
ID** 
KBAS 
NSTR 
VARY 
VARY 
VARY 
KARY 
KARY 


2 177 0.0 
2 178 0.0 
2 179 0.0 
2 180 0.0 
2 181 
2 182 
2 183 1.0 

2 184 

3 


NET C02 ADDITION (LB/HR) m n/HR ) 

OFT trace CONTAMINANTS ADDITION ( LB/HR ) 

. ... 

MISSION TIME (SEC) - 

MISSION TIME CURVE MULTIPLY FACTOR 

sand ABSORB CYCLE TIME r prjntoFF 

*1* PRINTOFF FREQUENCY, TIME STEPS FEK PKxn 

SPLIT.- MAIN. CABIN HUMIDITY CONTROL HX OR TO SAWO-I _ 


3 

3 

3 

4 
4 
4 
4 
4 

4 

5 
5 
5 

5 

6 
6 
6 
6 
6 
6 
6 
6 


10 


65 0.04866 SPLIT RATIO - GP0LY1 CALC 

MAIN CABIN SENSIBLE HX FAN . ~ 

1 CABIN FAN^ VOLUMETRIC FLOW R^E ^CFM) 

FAN ON/OFF SWITCH ( 1.0=ON,0.0-OFF ) 

fan HEAT ADO IT ION (WATTS) 

APlirTOl e UFAT R 


CONTINUATION 
OF K AND V 
COMPONENT 
DATA 


ADDITION 

SENSIBLE HEAT REMOVAL ONLY 


76 2100. 

84 1.0 _ 

MAIN CABIN SENSIBLE HX - — -601 

4 * | COUNTERF LOW » PRI=G AS , SEC= LIQ ,._SS CALC ... 

^ — — — H y rdUNTERFLOW UA (BTU/hR-F) 

‘ 6 cOOLiNG FLUID BOUNDARY COND - CABIN SENSIBLE HX 


49 


COOLANT FLOW (LB/HR) -GPOLYl CALC 

COOLANT TEMP IF) 


UUU Lnn I lU.r ' ■ ' 


ID** 

7 

KBAS 

7 

NSTR 

7 

VARY 

7 6 

ID** 

8 

KBAS 

8 

NSTR 

3 

VARY 

8 

VARY 

8 

VARY 

8 

ID** 

9 

_ KBAS 

9 

NSTR 

9 

VART 

9 

VARY 

9 

VARY 

- 9 

VARY 

9 

VARY 

9 

VARY 

9 

VARY” 

~ 9 

ID** 

10 

KBAS 

10 

NSTR 

10 

^ VARY 

10 

VARY 

10 

10** 

11 

KBAS 

11 

~ NSTR 

11 

ID** 

12 

KBAS 

12 


"'I - 1350 . " 

2 60.0 



MAIN CABIN CONDENSING HX -801 

4 31 2 _ --- PRI=GAS, SEC=LIQ. SS CALCS 

2 ° COUNTERFLOW HX UA (BTU/KR-F) _ GFOLY1 CALC 

66 COALING FLUID BOUNDARY COND - CABIN COND HX 

49 0 1 ___ 

— COOLANT F LOW ( LB/HR ) 

COOLANT TEMP (F) 

SEPARAT0S 

' 4 SPECIFY INDIVIDUAL SPLIT RATIOS 

SPLIT RATIO - COND VAFOR 

SPLIT RATIO - COND LIQUID ... 

SPLIT RATIO OXYGEN 

SPLIT RATIO - NITROGEN 
IpLIT RATIO - CARBON DIOXIDE 

SPLIT RATIO - TRACE CONTAMINANTS 

«. i - t-v c oftt i I F LOW SI 


1 950. * 

2 45.0 

3 50.0 

SPLIT 
10 7 

1 

66 0.03 

67 1.00 

68 0.03 

69 0.03 

70 0.03 

71 0.03 
“72 0.03 


10 


65 

67 


SPLIT RAT iU - 7 « r *«i" 

. 03 SPLIT RATIO - SPECIAL FL0U 31 

WATER SEPARATOR - CABIN COND HX 

1 " 9 2 CALC EXIT TMEP GIVEN Q - 

- - HEAT ADDED (BTU/HR) - GPOLYl CALC 

CONDENSATE REMOVED '^f.^^^^OR BYPASS 
GASMIX - COND HX WATER SEPARAiOR » W^EyEP ' 

^ 9 2 ... - 


COND HX FAN 
23 


11 


12 


13 


11 2 
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Table 8 (Continued) 
Input to ESCM 


NSTR 

VARY 

VARY 

VARY 

ID** 

KBAS 


12 1 

12 76 300. 

12 84 1*0 

12 91 140. 

13 GASMIX 
13 6 


NSTR 

13 



ID** 

14 

SPLIT - 

MAIN 

KBAS . 

14 

10 

- — - 

NSTR 

14 



VARY 

14 

65 0.0 


10** 

15 

GASMIX 

- STE 

KBAS 

15 .... 

__ 6 

— - 

NSTR 

15 

Z 


ID** 

16 

SAWD BED 31 

KBAS 

16 

73 0 

0 

NSTR 

16 


.... 

KARY 

16 

16 

0 

KARY 

16 

17 ‘ 

21 

KARY 

16 

18 

0 

KARY 

16 

19 

5 

KARY 

16 

20 

15 

VARY 

16 

60 70.0 


VARY 

16 

61 O.S 


VARY 

16 

63 70.0 


VARY 

16 

64 0.800 


VARY 

16 

65 1270. 


VARY 

16 

66 0.314 


VARY 

16 

67 32.70 

. 

VARY 

“16 

68 0.319 


VARY 

16 

69 10.0 


VARY 

16 

70 


VARY 

16 

71 


VARY 

16 

72 



CABIN FAN VOLUMETRIC FLOW RATE <CFM> 
fan on/off switch (1.0=ON,0.0-OFF) 
FAN HEAT ADDITION (WATTS) 

SAWO-II AND COND HX FAN . _ - 

12 2 5 

CABIN TO SAWD BED #1, SAWD^BED #2 


SPLIT RATIO - GPOLY1 CALC 

_ . — — — r>»p in ocn 


#1 INLET 


14 2 _ 


15 2 


-26 2 


.27 


16.. 


ALTCOM 


17 


VARY 
VARY 
VARY 
' VARY 
VARY 
VARY 
VARY 
" VARY 
VARY 


16 

16 

16 

16 

16 

16 

16 

16 

16 


73 

74 1.500 

75 0.400 

76 0.604 

77 70.00 

78 8.500 

79 7.400 

80 2.550 

81 0.005 


VARY 

16 

82 

VARY 

16 

83 

VARY 

16 

84 0.120 

VARY 

16 

35 180.0 

VARY 

16 

86 2230. 

VARY 

16 

87 240. 

VARY" 

”16 

88 

VARY 

16 

89 

VARY 

16 

90 

VARY 

16 

91 

" VARY 

16 

92 

VARY 

16 

93 0.470 

VARY 

16 

94 0.330 

VARY 

16 

95 100000. 

“ ID** 

17 

SPLIT - 

KBAS 

17 

10 

NSTR 

17 



SAWD 


absorb/d'esopb fug To = absorb, i = desorb ) 

COMPLEMENTARY BED'S COMPONENT NO 
DESORBING FLOW REVERSAL FLAG- 1-YES, 0 

numeer of bed segments — - — 

downstream component number 
ambient temperature (F) 

CAH-TO- AMBIENT KA/X ( BTU ^I F ’ , nEG _ F , 

EXTERIOR SURFACE TEMPERATU EXTERIOR > T BTU/HR-F ) 

DENSITY OF DRY RESIN BED, <P CF) ... - 

RESIN BED VOID VOLUME FRACTION . 



ISiSSsCf 

AVG STEAOY STATE BED TEMPERATURE, (DEG-F) 

TOTAL DRY RESIN BED WEIGHT (LBM) 

TOTAL CANISTER WEIGHT (LBM) — 

— TOTAL WATER IN BED (LBM) 

Sit BED H20 BE L0 ADINg\lB H20/LB BED) 



iSoM^lS.GECOEO, _ 

""' BED INLET DEW POINT, DEG-F 

III INLET°C02 P ?ARn AL E PRESSURE , m-M 

BED EXIT C02 PARTIAL PRESSURE, MM-HG _ 

- TOTAL BEO EXIT FLOW, CFM 

TOTAL WEIGHT OF FOAM IN BED. LBM 
coprTFIC HEAT OF FOAM BTU/LB-F 

TIME TO BEGIN DEBUG PR INTOFF .SECONDS. __ 

BE ! « TO CCHO HX OR C02 RECOVERY 2o g 

16 2 2 


CONTINUATION 

OF K AND V 

COMPONENT 

DATA 
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Table 8 (Continued) 
Input to ESCM 


VARY 

17 

65 0.0 

ID*# 

IB 

GASMIX - 

KBAS 

ia 

6 

NSTR 

la 


ID** 

19 

SAWD BE! 

KBAS 

19 

23 

NSTR 

19 

003 

VARY 

19 

2 70. 

VARY 

19 

3 14.7 

VARY 

19 

6 0.78 

VARY 

19 

10 20.55 

VARY 

19 

11 71.37 

VARY 

19 

12 0 . 30__ 

VARY 

19 

76 15. 

VARY 

19 

84 1.0 

VARY 

19 

91 50. 

ID** 

20 

GASMIX _ 

KBAS 

20 

6 

NSTR 

20 

2 

ID#* 

21 

SAWD BE 

KBAS 

21 

73 0 

NSTR 

21 


KARY 

21 

16 

KARY 

21 

17 

KARY 

21 

18 

KARY 

21 

19 

KARY 

21 

20 

VARY 

21 

60 70.0 

VARY 

21 

61 0.8 

VARY 

21 

63 70.0 

VARY 

21 

64 0.800 

VARY 

21 

65 1270. 

VARY 

21 

66 0.314 

VARY 

2 1~ 

67 32.70 

VARY 

21 

68 0.319 

VARY 

21 

69 10.0 

VARY 

21 

70 

VARY 

21 

71 

VARY 

21 

72 

VARY 

21 

73 

VARY 

21 

74 1.500 

'■"'VARY 

21 

75 0.400 

VARY 

21 

76 0.604 

VARY 

21 

77 70.00 

VARY 

21 

78 8.500 

VARY 

21 

79 7.400 

VARY 

21 

80 2.080 

VARY 

21 

61 0.300 

VARY 

21 

82 

VARY 

“21 

83 

VARY 

21 

84 0.120 

VARY 

21 

35 180.0 

VARY 

21 

86 2280. 

~ VARY 

21 

67 240. 

VARY 

21 

88 

VARY 

21 

89 

VARY 

21 

90 

■“ VARY 

21 

91 

VARY 

21 

92 

VARY 

21 

93 0.470 


SPLIT RATIO - GPOLY1 CALC 


17 2 


-22 2 


FAN 


IS 2 

INPUT <3 ONLY 

CABIN GAS MIXTURE: INITIAL TEMP (F) 

4 TOTAL PRESSURE J PSI A ) 

H20 VAPCR FLOW (LB/HR) 

OXYGEN FLOW < LB/HR ) 

NITROGEN FLOW (LB/HR) 

C02 FLOW (LB/HR) ... - 

CABIN FAN VOLUMETRIC FLOW RATE (CFM) 
FAN ON/OFF SWITCH ( 1 . 0=ON,0 . 0=OFF J 
FAN HEAT ADDITION (WATTS) 

_STEAM OR CABIN GAS TO SAWD BED #Z INLET 

-14 2 26 2 


30 2 


#2 SIMULATION 
20 2 


ALTCOM 


21 2 
22 2 


1 

16 

_0 

5 

20 


ABSORB/DESORB FLAG (0 = ABSORB » 1 = DESORB) 
COMPLEMENTARY BED’S COMPONENT NO. 

DESORBING FLOW REVERSAL FLAG: l=YES> p=NO 

NUMBER OF BED SEGMENTS 
DOWNSTREAM COMPONENT NUMBER 
AMBIENT TEMPERATURE (F) 

CAN-TO- AMBIENT KA/X ( BTU/HR-F ) _ , — 

EXTERIOR SURFACE TEMPERATURE, (DEG-F) 

CONDUCTANCE BETWEEN BED AND EXTERIOR > (BTU/HR-F) 

HEAT OF C02 REACTION WITH RESIN BED ( BTU/LB ) 

SPECIFIC HEAT OF DRY IR45 RESIN* (BTU/LB -F ) 

DENSITY OF DRY RESIN BED. ( FCF ) 

RESIN BED VOID VOLUME FRACTION 

TIME REQUIRED TO COMPLETE STEAM DESORB (SEC) 

GAS STREAM TOTAL SENSIBLE HEAT CHANGE ( BTU/HR ) _ 

GAS STREAM TOTAL LATENT HEAT CHANGE (BTU/HR) 

TOTAL C02 ABSORPTION* + )/DESORPTIQN( - ) RATE ( PPH ) 
TOTAL H20 EVAPORATION + )/CONDENSATION( - ) RATE ( PPH ) 
ABSORPTION/DESORPTION RATE ( PFH OF C02/LB-BED/PSI ) 
“AVG. STEADY STATE DESORPTION STEAM RATE (PPH) 

STEADY STATE C02 REMOVAL EFFICIENCY 

AVG. STEADY STATE BED TEMPERATURE, (DEG-F) 

TOTAL DRY RESIN BED WEIGHT (LBM) _ 

TOTAL CANISTER WEIGHT (LBM) 

TOTAL WATER IN BED (LBM) 

TOTAL C02 IN BED (LBM) 

TOTAL BED H20 LOADING (LB H20/LB BED) 

““ TOTAL BED C02 LOADING ( LB C02/LB BED) 

SPECIFIC HEAT OF CANISTER (BTU/LB-F ) 

EXIT TEMPERATURE TO EK’D DESORPTION, DEG-F 

TIME TO END A3S0RPTICN , SECONDS _ ... 

TIME LAG FOR OUTLET TEMPERATURE, SECONDS 
BED INLET DEW POINT, DEG-F 
BED EXIT DEW POINT, DEG-F 
BED INLET C02 PARTIAL FRBSSURE, MM-HG 
” bed EXIT C02 PARTIAL PRESSURE, MM-HG 
TOTAL BED EXIT FLOW, CFM 
TOTAL WEIGHT OF FOAM IN BED, LBM 
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Table 8 (Continued) 
Input to ESCM 


VARY 

21 

94 0.330 

VARY 

21 

95 100000. 

ID** 

22 

SPLIT - 

KBAS 

22 

10 

NSTR 

22 

„ 

VARY 

22 

65 1.0 

ID** 

23 

GASMIX 

KBAS 

23 

6 

NSTR 

23 


ID** 

24 

SPLIT - 

KBAS 

24 

10 

NSTR 

24 


VARY 

24 

65 0.0 

ID** 

25 

GASMIX 

KBAS 

25 

6 

NSTR 

25 


ID** 

26 

SPLITj 

KBAS 

26 

"10 

NSTR 

26 


VARY 

26 

65 1.0 

ID** 

27 

H20 SUI 

KBAS' 

~ 27 

30 3 

NSTR 

27 

1101 

VARY 

27 

1 2.55 

VARY 

27 

68 168.3 

VARY 

27 

“69 166.3 

VARY 

27 

70 70.0 

VARY 

27 

71 

VARY 

27 

72 30.0 

VARY 

27 

~ 98 

VARY 

27 

99 


SPECIFIC HEAT OF FOAM BTU/LB-F 

TIME TO BEGIN DEBUG PRINTOFF, SECONDS. 

' _ nrrnUCOY 


21 2 


23 


SPLIT RATIO - GPOLY1 CALC 


-22 2 


17 2 


24 


SPLIT"- SAKO EXIT GAS TO PREHEAT OR TO C02 ACCUMULATOR 


23 2 

SPLIT RATIO - GPOLYl„CALC 

"“STEAM GENERATOR AND PR EH EAT ^STEAMS 
29 2 


15 


24 2 


26 


25 2 


SPLIT RATIO - GPOLY1 CALC 


15 


25 


2 

2 


VARY 
ID#* 
~KBA$ 
NSTR 
VARY 
VARY 
ID** 
KBAS 
NSTR 
VARY 
'VARY 
VARY 
ID** 
KBAS 
NSTR 
VARY 
ID** 
KBAS 
K5TR 
ID** 
KBAS 
NSTR 
“ KART 
VARY 
VARY 
VARY 
"VARY 
TABL 
TITL 


COMPUTE OUTLET FLOW IN GPOLY1 
INITIAL H20 FLOW (LB/HR) 

H20 TANK MAX CAPACITY J LB) 

H20 TANK INITIAL FILL t LB ) 

H20 TEMPERATURE (F) 

TANK VOLUME (FT 3) 

TANK PRESSURE (PSIA) _ — — — 

WATER USEO FROM TANK (LB) 

WATER ADDED TO TANK (LB) 

H20 FLOW (PPH) 


27 100 ^q 5 suPPLY PUMP FOR STEAM GENERATOR 


28 
28 " 
28 
28 
28 
29 
29 
29 
29 
29 

29 

30 
30 
30 

30 

31 
31 

31 

32 
32 
32 
32 
32 
32 
32 
32 

1 
1 


22 
0002 
79 1.0 
85 10.0 


27 0 


INPUT PUMP HEAT ADDITION 
PUMP ON/OFF FLAG (1.0=ON) 

1Q o HEAT ADDITION (WATTS) _ _ 

STEAM GENERATOR FOR SAWD BED DESORPTION 

17 28 2 CALC ELEC ENERGY REQUIRED 

,, : 0 0 DESIRED DEGREES OF SUPERHEAT (F) 

TEMP OF SATURATED STEAM (F) 

ELEC ENERGY REQUIRED (WATTS) 


29 


25 


“67 

68 


ELEC cNtKw i Kcuiui^i-u 

cpLJT - DESORBED GAS TO CABIN OR C02 ACCUMULATOR. 

Jo 24 2 __ _ - - - 

SPLIT RATIO - GPOLY1 CALC 

GASMIX - CABIN AIR AND SAWD AIR TO HUMD CONTROL HX 

6 32 r 19 z — — 


32 


C02 ACCUMULATOR TANK 


30 
01000 

16 

70 70. 

71 2.0 

72 30.0 
60 0.4644 
1 2 


30 2 

co 2 delivery;' o = ovbd; 

C02 TEMPERATURE (F) 

TANK VOLUME (FT3) 

C02 PRESSURE (PSIA) 

MASS OF C02 IN TANK TLB) 
16 0 LIN LIN 


31 


1 r C02 REDUCTION. 


CREWMAN METABOLIC RATE ( BTU/HR ) VS MISSION TIME (SEC) 


END OF 
K AND V 
COMPONENT 
DATA 

4 - 

START OF 
TABLE DATA 
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Table 8 (Continued) 
Input to ESCM 


VALU 

1 

10 

21 

.0 

2700. 

3600. 

7200. 

.18Q0E+05 

. 1980E+Q5 

VALU 

1 

11 

20 

450.0 

450.0 

520.0 

650.0 

650.0 

450.0 

VALU 

1 

12 

21 

. 2340E+05 

. 2520E + 05 

. 3960E+05 

. 4140E+05 

.4320E+05 

.466CE+05 

VALU 

1 

13 

2D 

450.0 

650.0 

650.0 

520.0 

450.0 

450.0 

VALU 

1 

14 

21 

.5040E+05 

. 5400E + 05 

.5580E+05 

. 8640E+05 


. . 

VALU ~ 

1 

15 

2D 

390.0 

520.0 

300.0 

300.0 



TABL 

10 

1 


2 

3 

LIN STP 



TITL 

10 

2 

PRESSURE CONTROL SYSTEM OPERATION 




TITL 

10 

3 

PRESSURE MODE 

VS MISSION 

TIME (SEC) 





TITL 10 4 =0 14.5 PSI CONTROLLER USED 


TITL 

10 

5 


= 1 8.0 

PSI CONTROLLER USED 




VALU 

10 

10 

21 

-8400. 

15300. 

427864. 




VALU 

10 

11 

2D 

0.0 

0.0 

0.0 





. 

TABL 

20 

1 


2 

6 

LIN LIN 




TITL 

20 

2 

PRESSURE CONTROL SUBSYSTEM REGULATOR FLOWRATE 

CURVE 


TITL 

20 

3 

02 

OPENING 

FLOWS (FROM 

JULY 1980 PCS TESTING 

AT 

JSC) 


TITL 

20 

4 

02 

OPENING 

FLOWS VS TOTAL PRESSURE (PSIA) 




VALU 

20 

10 

21 

0.0 

14.510 

14.565 14.666 


14.700 

100.0 

VALU 

20 

11 

2D 

10.0 

10.0 

0.85 0.2 


0.0 

0.0 

TABL 

21 

1 


2 

6 

LIN LIN 




TITL 

21 

2 

PRESSURE CONTROL SUBSYSTEM REGULATOR FLOWRATE 

CURVE 


TITL 

21 

3 

02 

CLOSING 

FLOWS (FROM 

JULY 1980 PCS TESTING 

AT 

JSC) 


TITL 

21 

3 

02 

CLOSING 

FLOWS VS TOTAL PRESSURE (PSIA) 




VALU 

21 

10 

21 

0.0 

14.605 

14.672 14.745 


14.819 

100.0 

VALU 

21 

11 

2D 

10.0 

10.0 

0.6 0.25 


0.0 

0.0 

TABL 

22 

1 


2 

8 

LIN LIN 




TITL 

22 

2 

PRESSURE CONTROL SUBSYSTEM REGULATOR FLOWRATE 

UJ 

i 

u 


TITL 

22 

3 

N2 

OPENING 

FLOWS (FROM 

JULY 1930 PCS TESTING 

AT 

JSC) 


TITL 

22 

4 

N2 

OPENING 

FLOWS VS TOTAL PRESSURE (PSIA) 




VALU 

"22 

To 

21 

0.0 

14.500 

14.510 14.563 


14.583 

14.640 

VALU 

22 

11 

2D 

67.0 

67.0 

25.0 7.0 


1.0 

0.5 

VALU 

22 

12 

21 

14.748 

100.0 





VALU 

22 

13 

2D 

0.0 

0.0 





TABL 

23 

1 


2 

8 

LIN LIN 




TITL 

23 

2 

PRESSURE CONTROL SUBSYSTEM REGULATOR FLOWRATE 

CURVE 


TITL 

23 

3 

N2 

CLOSING 

FLOWS (FROM 

JULY 1900 PCS TESTING 

AT 

JSC) 


TITL 

21 

3 

N2 

CLOSING 

FLOWS VS TOTAL PRESSURE (PSIA) 




VALU 

23 

10 

21 

0.0 

14.590 

14.625 14.658 


14.680 

14.743 

VALU 

23 

11 

2D 

67.0 

67.0 

11.0 7.0 


0.8 

0.5 

VALU 

23 

12 

21 

14.813 

100.0 





VALU 

23 

13 

2D 

0.0 

0.0 





” TABL 

30 

f 


2 

12 

0 LIN LIN 





TITL 30 2 RELATIVE HUMIDITY (DECIMAL) VS ABSORPTION CYCLE DURATION (MIN) 


TITL 

30 

3 

FAN 

CHARACTERISTIC ADJUSTED - HSD 

REPORT SVHSER 

8921 


VALU 

30 

10 

21 

.15 

.25 

.35 

.45 

.50 

.55 

VALU 

30 

11 

2D 

13.00 

28.14 ~ 

41.00 

*51.57 " 

56 .00 

59.85 

VALU 

30 

12 

21 

.60 

.70 

.80 

.85 

.95 

1.C0 

VALU 

30 

13 

2D 

63.14 

68.00 

70.57 

70.99 

71.01 

71.01 

TABL 

31 

1 


2 

7 0 

LIN 

LIN 



TITL 

"31 

2 

SAWO FAN FLOW 

C CFM) VS ELAPSED TIME (SEC) DURING PREHEAT 


TITL 

31 

3 

ESTIMATED VALUES 1-9-85 





VALU 

31 

10 

21 

.0 

60. 

120. 

180. 

240. 

300. 

VALU 

31 

11 

2D 

1.90 

1.90 

3.50 

4.90 

5.90 

6.70 

VALU 

31 

12 

21 

"360. 






VALU 

31 

13 

2D 

7.10 






TABL 

32 

1 


2 

10 0 

LIN 

LIN 




TITL 32 2 SAND FAN FLOW ( CFM ) VS. ELAPSED ABSORPTION CYCLE TIME (SEC) 

“VALU 32 1C 21 .0 6C . 120. 240. 360. 450. 

VALU 32 11 2D 0.10 5.00 7.90 11.00 12.80 13.95 

VALU 32 12 21 600. 720. 900. 7200. 
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Table 8 (Continued) 
Input to ESCM 


VALU 

32 

13 20 14 

.60 

15.00 

15.00 15.00 



TABL 

33 

1 2 


10 0 

LIN LIN 



TITL 

33 

2 ABSORPTION CYCLE TIME (SEC) VS. SAWD FAN FLOW (CFM) 


VALU 

33 

10 21 0. 

10 

5.00 

7.90 11.00 12.80 

13.95 

VALU 

33 

11 2D .0 


60. 

120. 240. 360. 

430. 

VALU 

33 

12 21 14 

.60 

15.0 

15.0 15.0 



VALU 

33 

13 20 600. 

720. 

900. 7200. 



TABL 

41 

1 2 


16 0 

LOG LIN 



TITL 

41 

2 MOLE 

FRACTION 

I OF C02 IN 

BED SEGMENT VERSUS LOADING. 


VALU 

41 

10 21 

0.0 

0.01 

0.02 0.03 

0.04 

0.05 

VALU 

41 

11 20 

0000001 

0.000052 

0.00048 0.0019 

0.C050 

0.0105 

VALU 

41 

12 21 

0.06 

0.07 

0.075 0.0600 

0.0325 

0.085 

VALU 

41 

13 2D 

0.0225 

0.044 

0.068 0.1100 

0.15C0 

0.250 

VALU 

41 

14 21 

0.0885 

0.090 

0.094 1.0 



VALU 

41 

15 2D 

0.5100 

0.630 

1.000 1.0 



TABL 

51 

1 2 


12 0 

LIN LIN 



TITL 

51 

2 TEMPERATURE FACTOR ON CALCULATION OF EQUIL YC02. 



VALU 

51 

10 21 

55. 

77. 

90. 110. 

130. 

150. 

VALU 

51 

11 2D 

1.20 

1.00 

0.670 0.790 

0.570 

0.4400 

VALU 

51 

12 21 

170. 

190. 

210.0 250. 

270. 

400. 

VALU 

51 

13 2D 

0.2750 

0.1550 

0.060 0.0300 

0.0150 

0.0001 

TA3L 

61 

1 2 


12 0 

LIN LIN 



TITL 

61 

2 WATER 

LOADING 

FACTOR ON 

CALCULATION OF EQUIL YC02. 



VALU 

61 

10 21 

0.0 

0.09 

0.105 0.125 

0.157 

0.175 

VALU 

61 

11 2D 

0.551 

0.550 

0.571 0.745 

1.034 

1.200 

VALU 

61 

12 21 

0.200 

0.225 

0.250 0.275 

0.3G0 

1.000 

VALU 

61 

13 2D 

1.360 

1.410 

1.429 1.440 

1.450 

1.450 

TABL 

71 

1 2 


13 0 

LIN LIN 



TITL 

71 

2 INFLUENCE OF 

RELATIVE HUMIDITY ON EQUILIBRIUM WATER LOADING 

VALU 

71 

10 21 

-10.00 

0.000 

0.010 0.035 

0.060 

0.033 

VALU 

71 

11 2D 

0.00 

0.000 

0.100 0.200 

0 . 300 

0.400 

VALU 

71 

12 21 

0.104 

0.145 

0.190 0.230 

0.275 

0.315 

VALU 

71 

13 2D 

0.500 

0.600 

0.700 0.800 

0.900 

1.000 

VALU ' 

71 

14 21 

1.000 




0.315 

VALU 

71 

15 2D 

1.000 




1.000 

PLOTO 





ESCM DEMONSTRATION MODEL 



PLOT1 


2 181 



MISSION TIME (SEC) 



PLOT! 


2 182 



MISSION TIME (MIN) 




PLOT2 


1 66 



CREWMAN SEN ( BTU/HR ) 



PLOT2 


1 67 



CREWMAN LAT (BTU/HR) 



PLOT2 


1 63 



02 USE RATE (LB/HR) 



PLOT 2 


1 69 



C02 GEN RATE (LB/HR) 



PLOT2 


1 82 



CREWMAN MET ( BTU/HR ) 



PLOT2 


2 89 



REL HUMIDITY (DEC) 



PLOT 2 


2 94 



02 PRESSURE (PSIA) 



PLOT 2 


2 95 




N2 PRESSURE CPSIA) ~~ ' 

— — 

- - - 

PLOT2 


2 98 



DEW POINT TEMP (F) 



PLOT 2 


2 100 



C02 PRESS (MM HG) 



PLOT2 


2 104 



DRY BULB TEMP (F) 



PL0T2 


2 106 



"TOTAL PRESS (PSIA) 



PLOT2 


2 165 



02 MAKEUP (LB/HR) 



PLOT2 


2 166 



N2 MAKEUP (LB/HR) 



PLOT2 


15 6 



SA-1 H20 IN (LB/HR) 



~ PL0T2 


16 6 





5A-1 H20 OUT f LB/KR ) ~ 

. — — . 

— 

PLOT2 


15 12 



SA-l C02 IN ( LB/HR ) 



PLOT2 


16 12 



SA-1 C02 OUT (LB/HR) 



PLOTS 


16 2 



SA-1 EXIT TEMP (F) 



“ PLOT2 


16 82 



SA-1 H20 LCAD( LB/LB) 



PLOT2 


16 83 



SA-1 C02 LOAD( LB/LB ) 



PLOT2 


20 6 



SA-2 H20 IN (LB/HR) 
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START OF 
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PL0T2 

21 

6 

PLOT2 

ZO 

12 

— PLOT2 

21 

12 

PLOT2 

21 

2 

PLOT2 

21 

82 

PLOT2 

~ 21 

33 

EHDC 



ENDR 




SA-2 H20 OUT (LB/HR) 
SA-2 C02 IN (LB/HR) 
SA-2 C02 OUT (LB/HR) 
SA-2 EXIT TEMP IF) 
SA-2 H20 LOAD( LB/LB) 
SA-2 C02 LOAD( LB/LB) 


END OF 
PLOT DATA 
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Figure 8 

Input Setup And Case Data Loadsheet 
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m i ? i Case Data Group (Continued) 

4 * 1,4, format 

Record 1 

2) Print selection for Settle model 
using G189A. Not used by ESCM. 

Input a zero in column 8. 

IPR T = 0 =print phase instantaneous and 
average data , 

jPRT=l=print only instantaneous 
IPRT=2=delete all diagram printouts, 
do not print average data. 

see Figure 8 for the setup of the TAPE, BAS.C, and CASE fines of 
input. . ,,- 

In the case group data the 

of time in seconds ,TIMEMX, for value of 15 seconds has been 

^fectTto'gTve the f^tesTru^ning tin* while still givtng the 
desired computational accuracy. 

In the PR0P1 NAMELIST, the special fluids array are: 

1) liquid -ter at ir F and o l atmosphere 

2) carbon dioxide gas at /u r 

The diluent is nitrogen and the condensable is water vapor. 


4 . 1 . 2. 2 


if and V Component_Data_Grou£ 

ail the data required 

The K and V component .J ata ^ OU J 0 Monent. These data include the 
to specify and describe each c0 J p ?" e ? t er format and the numer- 
connection* and option djt. f «t«d The for-r are 

referred^ To ’"ill '*“ *° " 

V array values. 

. roniiired for each component. 

Five different types °f input are ^M r the KBAS record, the 
These input are contained in The IDrec ’ ds# jhe ID record 
NSTR record, the KARV records, an the user and describes what 

specifies the component number give y # of the crew is 

tC P»r^„WV foM^^ Vd. gives the format for 

the entries required. 
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4. 1.2. 2 K and V Component Data Group (Continued) 

The second record is the KBAS record. This input provides the basic 
input to describe how the components are connected together via 
fluid flow paths and also provides the solution sequence (i.e., the 
order in which each component is solved). The format for input to 
the KBAS record is shown in Figure 10. The input in Table 8 for the 
KBAS records are for the model shown in Figure 4 and are formatted 
per Figure 10. 

The third record in the K and V component group data is the NSTR 
record. In this record are input the component specific instruc- 
tions or options selected that are required by the component sub- 
routine. See Figure 11 for the format. A description of the 
options are provided in Reference 3 for each subroutine. Also, 
NSTR (16) through NSTR (18) have general meanings applicable to all 
subroutines. These general meanings can be found in Chapter 6 of 
Reference 3. 


The fourth and fifth records to be provided are the KARY and VARY 
records. These input provide the initial values and constants for 
integer and floating point variables required for each component. 
The format for inputting these values is given in Figures 12 and 13 
respectively for the KARY and VARY records. 


All the input for the K and V component group data are constants 
and should not be changed by the user. The only input that may be 
changed without affecting the ESCM model and correlation are the 
following: 


Item 


Location 


1) Number of crewmen 

2) Sensible heat per crewman 

3) Latent heat per crewman 

4) Cabin heat load 

5) Cabin free volume 

6) Printoff frequency 

7) Cabin gas design temperature 

8) Time to begin debug printoff 

9) CO 2 delivery: 

0=0 o bed 
l=Cu 2 reduction 


K( 1 ,16) 

V ( 1 ,66) 

V ( 1 ,67 ) 

V(2 ,66) 

V(2,139) 

V(2,184) 

V(2,87) 

V( 16 ,95)/V( 21 ,95) 
K(32,16) 


Other parameters may of course be changed, but only by a user who 
understands the ESCM program and the implications of the change. 
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4. 1.2. 3 Table Data Group 


4.2 


The table data group contains all the tables or curves required for 
the ESCM program. These tables relate such variables as metabolic 
inaH a function of time. These table data are not to be altered 
by the user except for crewman metabolic rate (Btu/hr) vs mission 

time (sec). 

Of course other table data may be changed by the user but should 
ly e Te done so by a user who understands the program and the 
implications of the change. The format to input 
shown in Figure 14. 


table data is 


4. 1.2. 4 Plot Data Group 


The plot data group of input specifies those variables which are to 
be saved for later use by a plotting program. The variables 1 sted 
in Table 8 serve as a baseline group from which some can be deleted 
or others added as desired by the user. The format for each vari- 
able to be entered is shown in Figure 15. 

Installation of HS ESCM Program at L angley Research Center 


4.2.1 Introduction 


The Hamilton Standard (HS) version of the ESCM program is 
at HS on IBM 3080 computers, and several modifications are 
make it compatible with the Langley Research Center (LaRC ) 
computer. These modifications involve substituting equival 

system functions (for example, date and c JJ tr *V.PJ°poJTRA N 
the IBM calls and updating the program to the local FOR IRAN 
procedure for making these changes to the program as deli ve 
is outlined below. All programs and procedure files are 
the directory SEB>LFR.DIR>EMSIM.DIR. 


executed 
needed to 
Prime 850 
ent prime 
time) for 
77. The 
red by HS 
found in 


4.2.2 Program Delivery 


HS will provide a nine-track magnetic tape with the prog )ram source 
code and input data at least. In addition listings of the code and 
results of a sample simulation run will needed _to verify -the 
correctness of program execution. This tape* QIC , word 
description of its block size, density, code (probably EBCDIC ) , i word 
length^ and number of files should be taken to the PRIME analyst w o 
will dump the tape to a user directory. 
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4.2.3 


Modification and Checkout 

The five subprogram calls that must be satisfied in the ESCM program 
are: DATEJC , JOBTIM, CORE, RTIME, and ERASE. These subprograms are 

not delivered with the program in any form since they are , IBM -system 
dependent. ERASE (called from ARSS) clears the interactive screen 
before data is to be displayed. This equivalent function has been 
provided by ERASE .F 77 and assumes that the terminal is a Tektronix 
terminal. The library VLBTEK.L must also be loaded to satisfy entry 
points required by the LaRC subroutine. DATEJC (called from ECLST) 
provides the date on the format MM/DD/YY. The equivalent function 
has been provided by DATEJC. F77. These two subroutines should be 
added to the delivered source code for the ESCM program which is i 
the file ESCM.F77. 

The three remaining functions are called from subroutines as 
JOBTIM from ELAPSE, CORE from ENKODE (which has an entry DEKODE), 
and RTIME from PRTIME. In these cases the calling subroutines 
themselves (which are included in the delivered source code) are 
to be replaced by source code of the same name, i.e., tiAHbt. , 
ENKODE. F77, and PRTIME. 77. These subroutines require that the prime 
library VAPPLB be loaded. 

The procedure file needed to compile this modified source code is 
called ESCM.CMPL.CPL and is listed below. To execute it, simply do 
R ESCM.CMPL. 


/* REPLACE SUBROUTINES ELAPSE, PRTIME, AND 
/* ENKODE WITH THE LRC VERSIONS 
/* ADD SUBROUTINES DATEJC AND ERASE 


F 7 7 ESCM.F77 - 64V - SILENT - DEBUG -XREF - LOUT 

This will produce a compilation listing on file OUT with variable 
cross references. FORTRAN errors, if any are encountered, will 
appear in the listing except for warning (level 1) messages (which 
can be numerous). When all FORTRAN errors have been corrected, this 
procedure will produce a binary file ESCM .BIN which can be loaded by 
procedure ESCM. LOAD. CPL listed below. To execute it, simply do 
ESCM. LOAD. 


R 


/* THIS LOADS THE ESCM BINARY AND LIBRARIES 

SEG -LOAD 

LO ESCM.BIN 

LI VAPPLB 

LI VLBTEK 

LI 

MAP LM 0 

SA 

Q 
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4.2.3 Modification and Checkout (Continued) 

This will produce a full load map on file LM which can be reviewed 
for unsatisfied entry points. If the load is complete, this will 
produce an absolute run file called ESCM.SEG. The program can be 
executed simply by SEG.ESCM. If errors are found, the same file can 
be run using the prime debugger by DBG.ESCM. Once all errors are 
corrected and an appropriate check case is verified, the -DEBUG 
option should be removed from the F77 command line in procedure 
ESCM.CMPL.CPL and both procedures re-executed to provide a more 
efficient run file. 

4.2.4 File Usage 

The ESCM program requires about a dozen files, several of which are 
site dependent and will not be used (on purpose) at LaRC. A summary 
table of the files is given below. These are defined by OPEN 
statements in program MAIN. File numbers are usually identified by 
integer variables which are usually set by DATA statements. Values 
for the LaRC prime version of ESCM are reflected in Table 9. 

4.3 ESCM Output 

Three forms of output are available with the ESCM program - output 
to the screen, hardcopy output, and plots. 

Output to the screen is automatic and will occur at the printoff 
frequency specified by Vocation V(2,184). For example, if the 
printoff frequency specified is 8.0 time steps per printoff, then a 
printoff will be made every 120 seconds when the recommended time 
step of 15 seconds is used. Values are printed in tabular form on 
the screen as a function of time. The values printed to the screen 
are: 

1) Time of operation, minutes 

2) Cabin temperature, °F 

3) Cabin dew point temperature, °F 

4) Cabin relative humidity, % 

5) Cabin CO 2 partial pressure, mmHg 

6) SAWD bed #1 CO 2 loading, % of dry amine wgt 

8) SAWD bed #1 exit temperature, °F 

9) SAWD bed #2 CO 2 loading, % of dry amine wgt 

10) SAWD bed #2 H 2 O loading, % of dry amine wgt 

11) SAWD bed #2 exit temperature, °F 
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ESCM Output (Continued) 

The hardcopy output is shown in Appendix B. It includes. 

1) A mirror image of the input. 

2) A more readable version of the entire input. 

3) Schematic printouts of the requested printoff frequency. See 
Figure 6 or Appendix B. 

4) Final values in the K and V arrays for all components at the end 
of the case. 

5) A list of variables whose time varying data have been saved for 
use by a plotting program. 

Plot output must be generated t ^°^ t gj e U 5 ^g°the n MERIAM P program 

which "wif^i ^turn "generate a'tape for use on the CALCOHP plotter. 
t generate plots, follow the following steps: 


1) LG EPLOT CNTL S0(G15UARL) RON 

2) LGEF GPLOT DATA G15UARL 


Make changes to NAMELIST CASE as required. The definition 

Ct d SIS 

setup S for* 1 t he™i ndenpendent $ and 

respectively. A sample setup is shown in Appendix B under 
GPLOT .DATA. 


A sample set of plots 


3) SUB EPLOT 

Execution of these steps will generate plots, 
is shown in Appendix D. 

Error Messages 

Listing of error messages can be found ^n ^sCM^re 

S&S3* 
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Table 10 

NAMELIST CASE Variables To Make Plots 



FORTRAN 

Format 

FORTRAN 

Name 

Definition 

— 

I 

NCASE 

G189 case number (Output by G189). 


I 

ILRCD 

Length of plot tape records (Output by G189). 

— 

I 

NPPTS 

Number of tape time points written for the G189 case. 
(Output by G189). 


E 

TSTART 

Plot start time (Seconds). 


E 

TEND 

Plot end time (Seconds). 


I 

IPFREQ 

Plot frequency flag. (Default value 3 l) 

— 

I 

NUMCAS 

Number of G189 tape cases flag. If more than one G189 
case is on the input plot tape, set NUMCAS=1. If only 
one G189 case is on the plot tape, NUMCAS*0. (Refer to 
G189 Input Data Setup for further details.) 

— 

I 

IFPLOT 

Plot data table generation flag. Must be equal to 0 
for the first plotting case. For succeeding plotting 
cases, set IFPLOT *1 if new plot data cards are not to 
be input. 

— 

I 

IEND 

End of run flag. 




IEND a 0 If this is not the last plotting case to 

be executed. 




IEND*1 If this is the last plotting case to be 

executed . 

— 

I 

IPLOT 

Plot only option flag. 




IPL0T=0 Plot editing information will be input. 




IPL0T=1 Plot only. No plot editing information 
will be input. (IFPLOT, LOCREF ,NVC) 


I 

IPRINT 

Plot file print flag. 


IPR I NT =0 No plot file time point history will be 
output . 

I PR I NT =1 A plot file time point history (list of 
time points for which a plot data record 
has been written) will be output. 
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Table 11 

Independent Variable Input To Make Plots 


Card Card FORTRAN FORTRAN 

Name Column Format Name 


Defi nitions 


PLOT 1-4 

5 II 


13-16 14 


19-28 E10.0 
29-38 E10.0 
39-80 7A6 


NVC* 


LOCREF* 


XMIN 

XMAX 


Variable type code. NVC=1 for 
the independent variable. If no 
independent variable card is 
input, the default independent 
variable is mission time (hours). 

Tape reference location for the 
independent variable. [Default 
independent variable is time 
(hrs). Tape Location 1.] 

Minimum value of the independent 
variable to be plotted (optional 
input) . 

Maximum value of the independent 
variable to be plotted (optional 
i nput ) . 


XI NTIT (7) Abcissa title. This alphanumeric 
information will appear on the 
plot frame as the horizontal axis 
title. 


* Do not input if I PLOT =1 
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Table 12 

Dependent Variable Input To Make Plots 


Card 

Name 

Card 

Column 

FORTRAN 

Format 

FORTRAN 

Name 

Definitions 

PLOT 

1-4 

5 

11 

NVC* 

Variable type code. NVC*2 for 
the dependent variable(s). 


6 

11 

NC** 

Number of dependent variable 
curves to be plotted for a frame. 
Default value is 1 (1 £ NC <_ 6). 


13-16 

14 

locref** 

Tape reference location for the 
independent variable. A table 
of tape reference location vs. 
V-array location (component number, 
reference location) is printed by 
G189. (Do not input on this card 
if NC > 1. Refer to multivariable 
data cards.) 


19-28 

E10.0 

YMIN 

Minimum value of the dependent 
vari ablets) to be plotted (optional 
input) . 


29-38 

E10.0 

YMAX 

Maximum value of the dependent 
variable(s) to be plotted (optional 
input) . 


39-80 

7A6 

YTITLE(7) 

Ordinate title. This alphanumeric 
information will appear on the 
plot as the vertical axis title. 

* Note* When using multivariable plots, the data value range of the 
dependent variables to be plotted on the frame should not 
differ significantly. 

** Do 

not input 

if IPL0T.1 





Multivariable plot Data Input 

Card 

Name 

Card 

Column 

FORTRAN 

Format 

FORTRAN 

Name 

Definitions 

PLOT 

1-4 

5 

11 

NVC 

Variable type code. NVC*2 for 
dependent variable. 


13-16 

14 

LOCREF 

Tape reference location for the 
ith dependent variable on a frame 

(i < 6). 


39-58 

5A4 

DVLAB(I) 

Label for the ith dependent vari- 
able on a plot frame. The varia- 
ble label and corresponding -plot 
character information are printed 
at the bottom of the plot frame. 

These 

cards are 

input only if 

NC > 1 and 

IPL0T*0. Do not input if IPL0T*1. 
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Error Messages (Continued) 


Subroutine 

SUB2 

SUB 2 


Message 


Subroutine SUB2 called with a temperature below 55°F ^ 
by component # Temperature — 

Subroutine SUB? called with a temperature above 300°F 
by component # Temperature - - 

outsfde^h^range^f^^e c^oi^rel a ti^on^for^the^ temperature eHect^n 
bed CO 2 loading. 


Subroutine 

BALNCE 


Message 

Subroutine IR45 called by component # 

Heat balance did not converge 

iterations. The time is seconds. 

energies = • Temperature = • 


, segment 

after 

Sum of the 
Vapor mass = 


This message alerts the h user that 

“ance'"^ energy eguatK.n This i message . may 

T no\ reach, ng 

convergence is less than O.b r. 

Lastly, although not an error, the user ^^p^y^ng^th’s TlT" 
location gs'^eHher'sAUD' b<5Tnput. The information printed is as 
fol 1 ows : 

Time 
K 

ICOUNT 
TBED 
TEMP 
TDEW 
XLOAD 

INERTS 

INERT 

PTOT 
PIF 

PH20 


Transient time, seconds 
Segment number 
Iteration number 
Bed temperature, °F 
Gas temperature, °F 

Dew point temperature, F , 

equivalent C0 2 loading, fraction of dry amine 
weight in time step DTIME, lbm 
Total inerts less C0 2 entering segment in time 

steD DTIME, 1b-mol6S • . 

: Total inerts less C0 2 initially in segment void, 

lb-moles . _ . _ 

: Total pressure of gas in segment, psia 

= Total pressure of gas in segment that .^°“ ld t * X s?J v 
if all gas entering segment were allowed to stay, 

= Equilibrium partial pressure of water, psia 
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Error Messages (Continued) 


PC02 
AE (I ) 

R2(I) 

B2(I) 

RATE 

PC021 

MCR1 

MCR2 

MC021 

MC02 

MH201 

MH202 

BURP 

PH20IN 

WTMAIN 

PC02IN 

PUSTRM 

PDSTRM 

QSUM 

HBEGIN 

HEND 

HC02 

HH20 

HSTOR 


= Equilibrium partial pressure of CO 2 , psia 
= Mass of gases exiting segment in time step DTIMt, 

1 bm . , 

= Mass of gases resident in void volume at end or 

time step, lbm 

= Mass of CO 2 in resin, H 2 O in resin, and dry resin 
at end of time step, lbm 

= C0 2 absorption/desorption rate, lb of C0 2 per lb of 
dry amine per psi difference between actual and 
equilibrium C0 2 pressure. 

= Initial partial rpessure of CO 2 from subroutine 
SUB2, psia 

= Initial calculation of mass of CO 2 absorbed using 
PC021 , lbm 

= Second calculation of mass of CO 2 absorbed after 
being limited to the masses of C0 2 actually 
present, lbm 

= Mass of CO 2 in gas entering plus resident after 
subtraction of amount absorbed MCR2, lbm 
= Mass of CO 2 in gas entering plus resident after 
absorption and after checks for pressure collapse 
and minimum CO 2 available to satisfy PCO 2 , lbm 
= Mass of H 2 0 in gas entering plus resident after 
evaporation or condensation but before pressure 
collapse check, lbm 

= Mass of H 2 0 in gas entering plus resident after 
evaporation or condensation and after pressure 
collapse check, lbm 

= Flag to indicate if bed temperature is at satura- 
tion temperature. T = true, F = false 
= Partial pressure of water vapor entering segment, 
psia 

= Molecular weight of gas mixture entering segment, 

1 b/1 b-mole 

= Partial pressure of C0 2 in gas entering segment, 
psia 

= Pressure upstream of segment, psia 
= Pressure downstream of segment, psia 
= Sum of all energies. When balanced QSUM = 0.0, Btu 
= Sum of energies in existing gas, entering gas, the 
bed, and tank at the beginning of the time step, 

Btu . . 

= Sum of energies in existing gas, exiting gas, tne 
bed, and tank at the end of the time step, Btu 
= Energy change due to absorption or desorption of 

CO 2 , Btu . 

= Energy change due to evaporation or condensation or 

water, Btu 

= Energy change due to the change in stored energy 
within the segment, Btu 
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4.4 Error Messages (Continued) 

MFG = Mass of steam condensed if positive, otherwise mass 

of steam evaporated. 

MCR = Mass of CO 2 absorbed if positive, otherwise mass of 

CO 2 desorbed. 

AI(I) = Mass of gases entering segment in time step 

DTIME , lbm 
I =l=total 

I=2=total non -condensable 
I=3=water vapor 
I =4=oxygen 
I=5=nitrogen 

r 1 ( I ) = Mass of gases resident in void volume at start, lbm 

B 1 ( I ) = Initial masses of (1) CO 2 in resin (2) H 2 O in resin 

(3) dry resin, lbm 

AIN(l) = Total mass of gas entering segment in time DTIME, 

lbm 

A I N ( 5 ) = Total non-condensables entering segment in time 

DTIME, lbm 

AIN(6) = Water vapor entering segment in time, DTIME, lbm 

AIN(9) = Molecular weight of non-condensable gases entering 

segment, lbm/lb-mole 

AIN( 12) = Carbon dioxide gas entering segment in time DTIME, 

1 bm 

This information is printed for each segment of the specified SAWD 
bed and each iteration in every time step following the user spec- 
ified time. 

5.0 FORTRAN NAME - ANALYSIS SYMBOL CROSS REFERENCE 

Table 13 is a cross reference list between the FORTRAN variable 
names used in the ESCM User's Manual and the analysis symbols used 
in the ESCM Model Description Document in Reference 4. Only those 
FORTRAN names that have a corresponding symbol in the Model Descrip- 
tion Document are listed. In some instances, the general symbol is 
presented for a specific FORTRAN name because that same symbol is 
used throughout the Model Description Document (i.e., m e = exiting 
mass flow, Ibm/hr). 

6.0 REFERENCES 

(1) Blakely, Robert L. and Rowell, Lawrence F.; "Environmental 
Control and Life Support System Analysis Tools for the Space 
Station Era"; SAE Technical Paper Series 840956; Fourteenth 
Intersociety Conference on Environmental Systems, San Diego, 
California; July 16-19, 1984. 

(2) Blakely, Robert "Contract NAS 1-17397, Development of an Emu- 
lation/Simulation Computer Model of a Space Station Environ- 
mental Control and Life Support System (ECLSS) Task 1 and Task 2 
Results", ESCM-EM-02, December 12, 1983. 
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Table 13 

FORTRAN Name - Analysis Symbol Cross Reference List 


FORTRAN 

Analysis 

Description 

Name 

Symbol 

CFM 

cfm 

Flow through fan, cfm 

DTIME 

At 

Time step increment, seconds 

KK( 1,16) 

N 

Number of people in crew 

R ( 1 ) 

m T 

Total flow exiting a component, lbm/hr 

R(2) 

T 

Temperature of fluid exiting a component, °F 

R(3) 

Pi 

Pressure at component inlet, psia 

R ( 4 ) 

Pe 

Pressure at component exit, psia 

R(5) 

m nc +m C02 

Total non-condensable flow, lbm/hr 

R(6) 

m v 

Vapor flow, lbm/hr 

R(7) 

m l 

Entrained liquid flow, lbm/hr 

R(8) 

C P 

Specific heat of non-condensables, Btu/lbm-°F 

R(9) 

M w 

Molecular weight of non-condensables, lb/lb-mole 

R ( io ) 

m 02 

Flow of oxygen, lbm/hr 

R ( 1 1 ) 

m N2 

Flow of nitrogen, lbm/hr 

R(12) 

m C02 

Flow of carbon dioxide, lbm/hr 

RHAVG 

RH 

Average relative humidity (decimal form) in cabin 
during absorption cycle 

V V ( 1 ,66) 

QSxN 

Crew total sensible heat generated, Btu/hr 

VV( 1 ,67) 

QLxN 

Crew total latent heat generated, Btu/hr 

VV(1,68) 

WO 2 

Crew oxygen usage rate, lbm/hr 

VV( 1 ,69) 

WC02 

Crew carbon dioxide generation rate, lbm/hr 

VV( 1 ,70) 

WH20 

Crew water vapor generation rate, lbm/hr 

V V ( 1 ,82) 

QT 

Crew total metabolic rate, Btu/hr 

VV (2,4) 

Pt 

Cabin total pressure, psia 

VV(2,87) 

^set 

Main cabin gas design or setpoint temperature, °F 

VV( 2,94) 

p 02 

Cabin partial pressure of oxygen, psia 

VV(2,165) 

mo2 

Flow of oxygen into cabin, lbm/hr 

VV(2 ,166) 

m N2 

Flow of nitrogen into cabin, lbm/hr 
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Table 13 (Continued) 


FORTRAN 
Name 

V V( 7 ,1) 

VV(7 ,66) 

VV(7,72) 

VV( 10,67) 

VV(10,65) 

VV(19,1) 

VV(32,71) 

VV(x, 65) 

The following apply to either SAWD bed: 

Specific heat of dry resin bed, Btu/lbm-°F 
Density of dry resin bed, Ibm/ft 

Basepoint absorption desorption rate (pph of C0 2 /lb 
of dry resin/psi of CO 2 partial pressure difference) 

Total dry resin weight, lbm 

Total thermal mass of tank or canister around SAWD 
resin, lbm 

Total bed water loading 
Total bed carbon dioxide loading 
Specific heat of canister, Btu/lbm-°F 
Outlet temperature sensor time constant, sec 
Total foam in bed, lbm 
Specific heat of foam, Btu/lbm-°F 


R(66) 

Cr 

R(67) 

'r 

R(74) 

r o 

R(78) 

M r 

R ( 79) 

Mt 

R(82) 

L H20 

R (83) 

l C02 

R(84) 

c t 

R(87) 

T 

R(93) 

M f 

R (94) 

Cf 


FORTRAN Name - Analysis Symbol Cross Reference List 


Analysis 

Symbol 


m a 

UA 

(Cp)eff 

mi 

Q 

(mj)e 

Vo 

SR 


Description 

Mass flow of air through humidity control heat 
exchanger, lbm/hr 

Calculated overall heat transfer coefficient of 
humidity control heat exchanger, Btu/hr- F 

Effective specific heat of gas How into humidity 

control heat exchanger, Btu/lbm-°F 

Condensate removal rate in water separator, lbm/hr 

Power added to water separator, Btu/hr 

Total gas flow exiting SAWD fan, lbm/hr 

CO 2 accumulator tank volume, ft 

Split ratio for splitter "x" 
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REFERENCES (Continued) 

( 3 ) "G189A General ied Envi ronmental /Thermal Control and Life 

Systems Computer Program Manual"; McDonnell Douglas Corporate 
MDAC-G2444; September, 1971. 

f4l Yanosv J "Model Description Document for a Computer Program of 

Report SVHSER 9504 for National Aeronautics and Space Admin 
tration Langley Research Center; NASA CR-181737, September 1988. 
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Appendix A 


Job Control Lists For IBM 


• ESCMCL 


• ESCMCN 
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ESCMCL 

JOB CONTROL LIST 


CONTROL PROMPT NOMSG 
Nl: DELETE PLOTL.DATA 
DELETE 01) TX .DATA 
FREE FILE( FT01F001 ) 

FREE FILE(FT02F001 ) 

FREE FILE(FT03F001) 

FREE FILE( FT04F001 ) 

FREE FILE(FT05F001 ) 

FREE FILE(FTOGFOOl) 

FREE FILE( FT10F001 ) 

FREE FILE<FT11F001 ) 

FREE FILE(FT16F001 ) 

FREE FILE(FT20F001 ) 

FREE ATTRL IST( A 1 B2 C3 D4) 
ATTRIB A1 BLKSIZE< 844 ) RECFM( V 
ATTRIB B2 BLKSIZE<5000) RECFM<V 
ATTRIB C3 BLKSIZEt 8404 ) RECFM( V 
ATTRIB D4 BLKSIZE<3990) RECFM(F 
DELETE DD1 
DELETE DD2 
DELETE DD3 
DELETE DD4 


B) LRECL<844) 

B S) LRECL(5000) 
B S) LRECL(8404> 
B A) LRECL< 133) 


DELETE DDA 

ALLO^DSCDDl) F(FTOIFOOI) NEW SPACE (844 ,200) BLOCK (844) ♦ 
USING(Al) CATALOG 

ALLOC DS( DD2 ) F(FT02F001) NEW SPACE(200,200) BLOCK(SOOO) + 
USING(B2) CATALOG 

ALLOC DS( DD3 ) F(FT03F001) NEW SPACE( 200,200) BL0CK(8404) ♦ 

USING( C3) CATALOG 

ALLOC DS(DD4) F(FT04F001) NEW SPACE (200, 200) BL0CK(840S) ♦ 
USING(C3) CATALOG 

££ S!SSSiS?i«S!SS!? , i» »««<*..*., . 

USING(D4) CATALOG 

ALLOC DS(DDA) F(FTIOFOOI) NEW SPACE(200,200) BL0CK(8408) ♦ 
USING(C3) CATALOG 

ALLOC DS(DDB) F(FTllFOOl) NEW SPACE(200,200> BL0CK(840o) ♦ 

ALLOC DS(PLOTl!dATA) F(FT16F001) NEW SPACE(200,200) BL0CK<8408) * 
UNIT(TSOWRKB) USING(C3) CATALOG 
ALLOC DS( * ) F ( FT20F001 ) 

CALL ’ ENG. G15 .LM( G1S9L2C) * 

END 





teSSologies 

M^R510[Lir®M 


Appendix B 


Files For Plotting 


• EPLOT 


• GPLOT 
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eplot.cntl 


PAGE 1 OF 1 


//TS0G15PD 

// 

// 

//* 

/♦JOBPARM 


J0 B <X.8888,99.08.,.,01..60.C189,1S603,0100,E>. 

>J TANOSTT ,N0TIFX-TS0G15P, 

CUSS-E.MSCLEVEL-(l,l),MSCCLASS-8 


LINES- 15 

/♦ROUTE PRINT LOCAL 
//GO EXEC PGM-GPLOT, REGION- 1024K 

//SIEPLIB OD DSNAME-ENG ' ^ ^ ' Ei^pgLETE > SPACE- ( 8408 , (209.290) ) . 

" FI02F 1S^ 

;7 T04F0 \cBMBEC F M-VBS.LKCL-8404 N0 W 

" r C n i;1sfG;Gf?LoT 0 i™!oilT-TSO*RKB.9ISP-(9L9.KEEP>, 

//ET16P001 99 , L RECL-8404,BLKSIZE-8408>. 

j / SPACE- (8408, (200,200) ,RLSE> 

//PETES EXEC MPRES, .MACII-B p . , olo , DELETE ) 

""£;££ »9 DSN“TS0G15P * TRES . DATA , 9ISP" ( 0L9 . DELETE ) 

//* 
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This appendix contains an ahnid 9 ed ootput fro* the fo,1»1n, sa m p,e pr*l- 
The principal input parameters are. 

(1) Number of men =3 

(2) Cabin volume - 8,0UU rt 

(3) Cabin temperature = 70 F 

(4) Equipment heat load - 17065 Btu/hr 

(5) Cabin leakage = 0.08333 lbm/ r t ati ve 

6 A variable metabolic load representative 

of a work, eat, sleep schedule 
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The following plots in this Appendix were obtained from the data generated by 
the sample problem of Appendix C. The Hamilton Standard MERIAM program was 
used to generate a tape for use on the CALCOMP plotter. The procedure des- 
cribed in Section 4.3 was used with the plotting files shown in Appendix B. 
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